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Preface 



Sustainable Development is defined as meeting the needs of the current generation 
without compromising the ability of future generations to meet their own needs. 
Nowadays, Sustainable Development has been the subject of wide-ranging dis- 
cussions and debates within government, non-government and academic circles, 
being a major focus of national and international economic, social and environ- 
mental agendas so that a good quality of life can be enjoyed by current and future 
generations. However, Sustainable Development is threatened by industrial pol- 
lution emissions which cause environmental problems such as deforestation, soil 
erosion, acid rain, toxicity, water pollution, depletion of fisheries and the expan- 
sion of desert areas. The problem of climate change, which occurs due to the 
increasing concentrations of greenhouse gasses in the atmosphere, has caused great 
concern at all levels, from the general public to national governments and 
international agencies. Renewable energy sources can be an important remedy to 
many environmental problems that the world faces today. Some new governmental 
policies have been adopted to encourage the introduction of energy efficiency 
measures, the technical changes and the renewable and sustainable energy. The 
energy planning scenario has completely changed over the las two decades, from 
an almost exclusively concern with cost minimization of supply-side options to the 
need of explicitly multiple and potentially conflicting objectives. In addition, 
different and numerous groups of actors, such as institutions and administration 
authorities, potential investors and environmental groups, get involved in the 
process. This complex environment indicates the multi-criteria character of the 
problem. The traditional single criteria decision-making approach is no longer able 
to handle these problems. The policy formulation for fossil fuels energy substi- 
tution by renewable energies must be addressed in a multi-criteria context. Multi- 
criteria decision-making is a decision support approach that is suitable for 
addressing complex problems featuring high uncertatinty, conflicting objectives, 
different forms of data and information and multi interests and perspectives. The 
multiplicity of criteria and the complexity of energy planning and energy projects 
makes multi-criteria analysis a valuable tool in the decision-making process. It is 
the aim of this book to show the use of the multi-criteria decision-making methods 
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to provide valuable assistance in reaching equitable and acceptable solutions in the 
selection of renewable energy projects. This book is aimed at those actors who 
bring along different criteria and point of view, have interests at stake and 
therefore, influence the decision-making process. 

Cantabria, Spain, November 2011 Jose Ramon San Cristobal Mateo 
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Chapter 1 

The Renewable Energy Industry and the Need 
for a Multi-Criteria Analysis 



Energy systems play an important role in the economic and social development of a 
country. The selection among different alternatives for the expansion of renewable 
energy sources, where the environmental dimension is added to the usual economic, 
technical and social criteria demands the utilization of powerful decision-aid tech- 
niques. The traditional single decision-making approach is no longer able to handle 
these problems. The multiplicity of criteria, and the involvement of different actors 
in the decision procedure, makes multi-criteria analysis a valuable tool in the policy 
formulation for fossil fuel energy substitution by renewable energy sources. 



1.1 Introduction 

Renewable energy (RE) is the energy generated from natural resources such as sun- 
light, wind, rain, biomass, tides and geothermal heat. The exploitation of RE sources 
has gained a vast interest during the last years. Global environmental problems, 
including the greenhouse effect and ozone holes, which have a negative impact on 
the health and living quality of humans, are now identified as significant problems. 
It is now widely recognized that the largest source of atmospheric pollution stems 
from fossil fuel combustion, upon which current energy production and use patterns 
throughout the world rely heavily [1]. The energy demand of fossil fuel has major 
consequences on the environment. The emission of toxic chemical pollutants, green- 
house gasses, and other pollutants cause climate change and environmental pollution 
of air, land and water. In addition, the growing rise in the cost of raw materials driven 
up by oil prices, which are strongly affected by political conditions in the world, can 
have a negative effect on the world economies. In this context, the raising aware- 
ness of environmental issues, the precarious nature of dependency on fossil fuel 
imports, the vulnerability of many states as regard energy sector and the advent of 
RE alternatives, have forced many countries, specially the developed ones, to use RE 
sources. 
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1 The Renewable Energy Industry 



Clean, domestic and renewable energy is commonly accepted as the key for future 
life [2]. Setting goals to ensure that the environment is safeguarded and the effects of 
climate change are mitigated while seeking long-lasting Sustainable Development 
is one of the most important issues faced by the world today [3]. The Sustainable 
Development has been the subject of wide-ranging discussion and debate within gov- 
ernment, non- government and academic cycles, being a major focus of national and 
international economic, social and environmental agendas. Sustainable Development 
is defined as meeting the needs of the current generation without compromising the 
ability of future generations to meet their own needs [4]. Energy systems play an 
important role in the economic and social development of a country and the living 
of people. A balance of energy production and consumption which has no, or mini- 
mal, negative impact on the environment, but gives the opportunity for a country to 
employ its social and economic activities can be seen as the final goal [5]. 

The use of renewable energy sources for the production of electric power brings 
huge benefits, both in terms of environmental protection as well as savings in non- 
renewable resources [3]. Renewable energy sources are environment-friendly com- 
pared to fossil fuels, with low levels or absence of both air pollutants and greenhouse 
gasses emissions, waste valuation potential and capable of replacing conventional 
sources in a variety of applications at competitive prices [6]. A more intensive use 
of renewable energy sources by means of improvements in energy performance and 
transport, conversion and conservation technologies, as well as development of new 
energy sources, can be an important remedy to many environmental problems that the 
world faces today. Some new governmental policies have been adopted to encour- 
age the introduction of energy efficiency measures, the technical changes and the 
renewable and sustainable energy. The Kyoto Protocol, the EU Renewable Directive 
2001/77/EC [7] and the European Biomass Action Plan [8] are examples of ambi- 
tious political goals fostering the development of conversion technologies based on 
Renewable energy sources. 



1.2 The Need for a Multi-Criteria Analysis 

Given that the energy sector and energy planning in particular affects the interests and 
resources of multiple actors, it is socially not acceptable to suggest (or even imple- 
ment) a policy alternative without taking into account the interests and preferences 
of the (multiple) affected factors [9]. Different and numerous groups of actors get 
involved in the process, each group brings along different criteria and point of view, 
which must be resolved within a framework of understanding and mutual compro- 
mise (concessions) [10]. The actors include those groups of individuals, institutions 
and administration authorities such as local authorities and communities, potential 
investors, academic institutions, environmental groups, governments, that through 
their priorities and evaluation systems, have interests at stake and directly or indi- 
rectly infiuence the decision-making process. 
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Table 1.1 Main criteria in the renewable energy industry [5] 



Criteria 


Description 


Technical criteria 


Efficiency 


It is the ratio of the output energy to the input energy 


Exergy efficiency 


Exergy is the energy that is available to be used 


Primary energy ratio (PER) 


PER is defined as the ratio of consumption primary 
energy to the users’ demand energy 


Safety 


Safety can be considered an evaluation criteria to show 
the effects of energy systems to society and people 


Reliability 


It may be defined as the capacity of a device or system 
to perform as designed 


Maturity 


It measures the degree of maturity of the technology 
showing how widespread the technology is at both 
national and international level 


Environmental criteria 


Particles emissions 


It involves mainly NO^, CO 2 , CO and SO 2 emissions 
which can lead to global warming and contribute to 
the greenhouse effect 


Non-methane volatile organic 


NMVOC is a generic term for a large variety of 


compounds (NMVOCs) 


chemically different compounds, like for example, 
benzene, ethanol, formaldehyde, cyclohexane or 
acetone 


Land use 


The environment and landscape are affected directly by 
the land occupied by energy systems. Quality of 
people’s life is affected by energy systems as it could 
have been used for creation of parks and recreation 
centers. Particularly energy supply systems with 
biomass and biofuels require a large amount of land 


Noise 


Noise pollution form energy systems is displeasing 
machine-created sound that disrupts the activity or 
balance of human and animal life 


Topography 


It may play an important role when selecting between 
several alternatives such as wind farm locations 



Another factor that has completely changed the energy planning scenario over 
the last two decades, from an almost exclusively concern with cost minimization 
of supply-side options to the need of explicitly multiple and potentially conflicting 
objectives is the increasingly complex social, economical, technological, and envi- 
ronmental scenario. Multiple criteria that affect the success of a renewable energy 
project must be analyzed and taking into account. The literature traditionally divides 
the used criteria into four aspects technical, environmental, economical, and social 
as shown in Tables 1.1 and 1.2 [5]. 

The traditional single criteria decision making approach, normally aimed at iden- 
tifying the most efficient option at a low cost, is no longer able to handle these prob- 
lems. From an interdisciplinary point of view, exploiting renewable energy sources 
(potential solutions, conrses of actions, decision alternatives, selecting the best 
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1 The Renewable Energy Industry 



Table 1.2 Main criteria in the renewable energy industry [5] 



Criteria 


Description 


Economic criteria 
Investment cost 


Purchase of mechanical equipment, technological 
installations, construction of roads and connections 
to the national grid, engineering services, drilling 
and other incidental construction work 


Operation and maintenance costs 


Operation cost such as employees’ wages, the funds 
spent for the energy, the products and services for the 
energy system operation and maintenance costs that 
aim to prolong energy systems life and avoid failures 


Fuel costs 


It refers to the fund spent for the provision of raw 
material necessary for energy supply system 
operation 


Electric cost 


It is the product cost of power plant, to evaluate its 
economic performance from the viewpoint of 
consumers 


Internal rate of return 


It is a commonly used investment criteria expressing the 
rate of return of a project 


Net present value 


It is defined as the total present value of a time series of 
cash flows 


Payback period 


Payback period refers to the period of time required for 
the return on an investment to ‘repay’ the sum of the 
original investment 


Service life 


It is the expected lifetime of an energy system 


Equivalent annual cost 


Equivalent annual cost (EAC). It is the cost per year of 
owning and operating an asset over its entire life 
span. It is calculated by dividing the NPV of an 
energy project by the present value of an annuity 
factor 


Social criteria 
Social acceptability 


It is a qualitative criteria expressing the overview of 
opinions related to the energy systems by the local 
population regarding the hypothesized realization of 
the projects under review from the consumer point of 
view 


Job creation 
Social benefits 


Number of jobs created during the project’s cycle life 
It expresses the social progress in the local regions that 
house the energy systems, estimating the total social 
and economic impact that may become perceptible 
by introducing an energy project 



alternative, identifying non-sustainable energy supply system, informing design- 
makers of the integrated performances of the alternatives and monitoring impacts 
on the social environment, concentrate investments in the most appropriate areas) 
can no longer be based on a single-dimensional axis of evaluation, such as cost or 
benefit. The policy formulation for fossil fuels energy substitution by RE must be 
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addressed in a multi-criteria context. The multiplicity of criteria and the complexity 
of energy planning and energy projects make multi-criteria analysis a valuable tool 
in the decision-making process. Decision-makers need to examine more than one 
objective in reaching the final goal set with respect to the satisfaction of constraints 
set by the environment. 
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Chapter 2 

Multi- Criteria Analysis 



Decision-making is the study of identifying and choosing alternatives to find the 
best solution based on different factors and considering the decision-makers’ expec- 
tations. Every decision is made within a decision environment, which is defined as 
the collection of information, alternatives, values and preferences available at the 
time when the decision must be made. The difficult point in decision-making is the 
multiplicity of the criteria set for judging the alternatives. The objectives are usu- 
ally conflicting and, in most of the cases, different groups of decision-makers are 
involved in the process. To facilitate this type of analysis, a family of tools referred to 
as Multi-criteria decision-making methods gained ground due to the need to have a 
formalized method to assist decision-making in situations involving multiple criteria. 

2.1 Introduction 

Multi-criteria decision-making methods is a branch of a general class of Opera- 
tions Research models that is suitable for addressing complex problems featuring 
high uncertainty, conflicting objectives, different forms of data and information, 
multi interests and perspectives, and the accounting for complex and evolving bio- 
physical and socio-economic systems [1]. This major class of methods is further 
divided into Multi-objective decision-making and Multi-attribute decision-making. 
These methodologies share the common characteristics of conflict among criteria, 
incommensurable units and difficulties in design/selection of alternatives. The main 
distinction between the two groups of methods is based on the number of alterna- 
tives under evaluation. Multi-attribute decision-making methods are designed for 
selecting discrete alternatives while Multi-objective decision-making methods are 
more adequate to deal with multi-objective planning problems, when a theoretically 
infinite number of continuous alternatives are defined by a set of constraints on a vec- 
tor of decision variables [2-4]. In Multi-objective decision-making (also known as 
multi objective programming or a vector optimization/maximization/minimization 
problem), the alternatives are not predetermined but instead a set of objective 
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2 Multi-Criteria Analysis 



functions are optimized subject to a set of constraints. The most satisfactory and 
efficient solution is sought. In this identified efficient solution it is not possible to 
improve the performance of any objective without degrading the performance of 
at least one other objective. In Multi-attribute decision-making a small number of 
alternatives are to be evaluated against a set of attributes which are often hard to 
quantify [5]. 

Following Tsoutos et al. [6] there are four starting reasons that justify the use 
of Multi-criteria decision-making methods: (i) It allows for investigation and inte- 
gration of the interests and objectives of multiple actors since the input of both 
quantitative and qualitative information from every actor is taken into account in 
form of criteria and weight factors; (ii) It deals with the complexity of the multi- 
actor setting hy providing output information that is easy to communicate to actors. 
The user-friendliness of the method lies on two aspects: the suggested criteria are 
estimated and given values that are consistent and comparable with the input data 
(as a measure of appropriateness); and the ‘simple’ format of the output of the method 
that makes the method’s results meaningful and directly applicable for the interested 
actors; (iii) It is well-known and applied method of alternatives’ assessment that 
also includes different versions of the method developed and researched for specific 
problems and/or specific contexts and (iv) It is a method that allows for objectivity 
and inclusiveness of different perceptions and interests of actor without being energy 
and cost intensive. 

These methods can provide solutions to increasing complex management prob- 
lems. They provide better understanding of inherent features of decision problem, 
promote the role of participants in decision-making processes, facilitate compromise 
and collective decisions and provide a good platform to understand the perception of 
models and analysts in a realistic scenario. The methods help to improve quality of 
decisions by making them more explicit, rational and efficient. Negotiating, quan- 
tifying and communicating the priorities are also facilitated with the use of these 
methods [7]. 

It should be noted that methods and results are not necessarily comparable. Every 
method has its restrictions, mostly due to model assumptions, which should be con- 
sidered when the method is used. Inconsistences might arise because [8]: (i) the 
choice problem formulations do not reflect the same preference structures, (ii) the 
ways in which preference information is processed vary between different methods 
and (iii) the methods interpret the criterion weights differently. 



2.2 Model Building 

Multi-criteria decision-making may be considered as a complex and dynamic process 
including one managerial level and one engineering level. The managerial level 
defines the goals, and chooses the final optimal alternative whereas the the engi- 
neering level defines the alternatives, points out the consequences of choosing any 
one of them from the stand point of various criteria and performs the multi-criteria 
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ranking of alternatives. The engineering level performs the optimization procedure. 
At the managerial level, public officials, called decision-makers have the power to 
accept or reject the solution proposed by the engineering level [ 9 ]. The decision- 
making process usually includes five main stages: defining the problem, generating 
alternatives and establishing criteria, criteria selection, criteria weighting, evaluation, 
selecting the appropriate multi-criteria method and finally ranking the alternatives. 
The main steps of Multi-criteria decision-making are the following: 

Step 1. Defining the problem, generating alternatives and establishing criteria 

A decision-making problem should start out by clearly defining the problem, discern- 
ing the alternatives, identifying the actors, the objectives and any points in conflict, 
together with the constraints, the degree of uncertainty and the key issues. After this, 
the problem can be framed indicating the evaluation criteria. 

Step 2. Assigning criteria weights 

The next steps include the assignment of criteria weights. These weights, that show 
the relative importance of criteria in the multi-criteria problem under consideration, 
can be determined by techniques such as Analytical Hierarchy Process and Sismos 
approach. 

Step 3. Construction of the evaluation matrix 

The phase in which the model is built constitutes a process from which the ‘essence’ 
of the problem is extracted from the complex picture drawn up so that the problem 
can be assessed adequately. At the end of this step, the MCDM problem can be 
expressed in matrix form as follows: 



Criteria Ci,C2,...,C„ 

Weights W\,W2, ■ ■ ■ ,Wn 

Alternatives 



■ Ar 




Xll X\2 ■■ ■ Xln 


A 2 




X2\ X22 ■■ ■ X2n 


Am _ 




Xm\ Xfyi2 * * • Xj^yi 



( 2 . 1 ) 



where x,y is the evaluation given to alternative / th with respect to criterion yth, iVj is 
the weight of criteria j, n is the number of criteria and m is the number of alternatives. 

Step 4. Selecting the appropriate method 

A multi-criteria method must be selected and applied to the problem under consid- 
eration in order to rank alternatives. The data and the degree of uncertainty are key 
factors for the decision-maker when selecting among several multi-criteria methods. 
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Step 5. Ranking the alternatives 

Finally, the alternatives’ ranking is ordered and the best ranked alternative is proposed 
as a solution. 
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Chapter 3 

AHP 



When, for a given problem, a decision-maker wishes to evaluate the performance of 
a number of alternative solutions and multiple objectives are important, it may be 
difficult to choose between alternatives. These alternatives can be evaluated in terms 
of a number of decision criteria and sometimes an alternative is superior in terms 
of one or some of the decision criteria, but inferior in terms of some other criteria. 
There are several reasons for adopting the AHP method as a multi criteria decision- 
making tool. For instance, qualitative criteria, such as political or social factors, need 
to be included in the decision-making process; a large quantity of criteria need to 
be considered or, indeed, when dealing with future technologies, characteristics of 
these technologies may not be well known. 



3.1 Introduction 

The Analytical Hierarchy Process (AHP) is a decision-making procedure which was 
originally developed by Saaty in the 1970s [1], and since then has been an effective 
tool In structuring and modeling situations involving multiple and sometimes con- 
flicting objectives. The objective of using AHP is to identify the preferred alternative 
and also determine a ranking of the alternatives when all the decision criteria are 
considered simultaneously. 

AHP is applied by breaking down an unstructured problem into components 
parts. Hierarchical orders are arranged by forming hierarchical structures with goal 
(objective) at the top of the hierarchy, criteria affecting the decision at the inter- 
mediate level, and decisions options (alternatives) at the bottom of the hierar- 
chy. The decision-maker is then guided through a series of pair-wise comparison 
judgments to express relative strength or intensity of impact of the elements in the 
hierarchy. These judgments, represented by assigning numerical values, are then 
synthesized in the use of eigen vectors to determine which variables have the highest 
priority. 
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3 AHP 



Table 3.1 Integer values and interpretation 



Oij Value 


Interpretation 


1 


Objectives / and j are of equal importance 


3 


Objective i is weakly more important than objective j 


5 


Experience and judgements indicate that objective i is strongly 
more important than objective) 


7 


Objective i is very strongly or demonstrably more important 
than objective j 


9 


Objective i is absolutely more important than objective j 


2, 4, 6, 8 


Intermediate values, for example, a value of 8 means that 

objective i is midway between strongly and absolutely more 
important than objective j 



AHP is an intuitive method for formulating and analyzing decisions which uses 
the following four steps in solving a problem: Structuring the decision problem into 
a hierarchical model, Obtaining the weights for each criteria, Finding the score of 
each alternative for each criteria and Obtaining and overall score for each alternative. 

Step 1. Structuring of the decision problem into a hierarchical model 

Each energy technology has many components and its use has social, environmental, 
political or economical implications. These are sometimes dependent on each other. 
The first step involves structuring of the decision into a hierarchical model built 
with support from specialists in applicable fields. This involves the decomposition 
of the problem into elements according to their common characteristics and forming 
a hierarchical model at different levels showing the relationships between the goal, 
criteria (subcriteria if necessary) and alternatives. 

Step 2. Obtainiug the weights for each criteria 

In the second step, once the relevant criteria and alternatives have been determined 
together with the relations between them, the elements of a particular level are com- 
pared pair-wise with respect to a specific element in the immediate upper level. A 
judgmental matrix, the pairwise comparison matrix (A), is formed using linguistic 
terms that include the verbal judgments of the decision-maker. The entry in row i 
and column j of A (called aij) indicates how much more important criterion i is than 
j with respect to the alternative. To transform the verbal judgments into numerical 
quantities, Saaty [1] suggested the use of integer values constituting a 1-9 scale as 
shown in Table 3.1. 

When constructing the pairwise comparison matrix, the following rules must be 
verified: 

1 . If Gij — a, then a ji = ^ . 

2. If criteria i is judged to be of equal relative importance as criteria j, then 
aij — Gji — 1, and an — 1 for all i. 
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Table 3.2 RI for different values of n 



n 


2 


3 


4 


5 


6 


7 


8 


9 


10 


11 


12 


RI 


0 


0.58 


0.90 


1. 12 


1.24 


1.32 


1.41 


1.45 


1.49 


1.51 


1.48 



3. If all the comparisons are perfectly consistent, then the relation a,* = atj — 
ajkii, j, k. 

To recover the vector W — [u>i , W 2 , ■ ■ ■ , u>n] from A, which indicates the weight 
that each criteria is given in the pairwise comparison matrix, the following two-step 
procedure is used: 

1 . For each of the As column divide each entry in column i of A by the sum of the 
entries in column i. This yields a new matrix, called Anorm (for normalized) in 
which the sum of the entries in each column is 1 . 

2. Estimate Wi as the average of the entries in row i of Anorm- 

Once we have the pairwise comparisons matrix it is necessary checking it for 
consistency. Slight inconsistencies are common and do not cause serious difficulties. 
We can use the following four-step procedure to check for the consistency on the 
decision-maker’s comparisons. 

1. Compute AW^ , where W denotes our estimate of the decision-maker’s weight. 

2. Find out the maximum Eigen value {kmax)'- 

1 ith entry in AVT^ 

-T (3.1) 

n /th entry m W‘ 

3. Compute the Consistency Index (Cl) as follows: 

Cl = ~ ” (3.2) 

M — 1 

The smaller the Cl, the smaller the deviation form the consistency is. If Cl is suf- 
ficiently small, the decision-maker’s comparisons are probably consistent enough 
to give useful estimates of the weights for their objective. For a perfectly consis- 
tent decision-maker, the ith entry in AW^ — n(ith entry of W^). This implies 
that a perfectly consistent decision-maker has Cl = 0. 

4. Compare the Consistency Index to the Random Index (RI) for the appropriate 
value of n, shown in Table 3.2. 

If ^ < 0.10, the degree of consistency is satisfactory, but if ^ > 0.10, serious 
inconsistencies may exist, and the AHP may not yield meaningful results. 

Step 3. Finding the score of each alternative for each criteria 

We now determine how well each alternative ‘satisfies’ or ‘scores’ for each objective. 
To determine these scores, a pairwise comparison matrix for each objective must be 
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Table 3.3 Criteria to evaluate 
a wind farm location 



Criteria 



Cl 

C2 

C3 

C4 

C5 



Topography 

Operation and maintenance cost 
Land use 
Infrastructure 
Investment cost 



constructed. In this matrix, the rows and columns are the decision-maker’s possible 
decisions. 

Step 4. Obtaining overall score for each alternative 

Synthesis is the process of weighting and combining priorities through the model 
that leads to the overall results. Through this process, we synthesize the objective 
weights with the scores of each alternative on each criteria to obtain an overall score 
of each alternative. For calculating the ‘composite weight’, the following equation 
is used: 

2] (Local priority of A 1 with respect to Cl) 
Final priority of alternative 1 = / 

X (Local priority of Ci with respect to goal) 

(3.3) 



3.2 Application 

In this section, the AHP is applied to the selection of the most convenient location 
for a wind farm project in the autonomous region of Cantabria, in the north of Spain. 
Four locations. La Bragui'a (Ai), Cilda (A 2 ), San Pedro (A 3 ) and Estacas (A 4 ) will 
be evaluated according to the criteria shown in Table 3.3. 

A typical AHP structure based on a three-level hierarchy is shown in Fig. 3.1 
where the top most level represents the main goal or focus of the problem. The 
intermediate levels correspond to the criteria and sub-criteria, while the lowest level 
contains the decision alternatives. 

In the selection of the RE alternatives, the following pairwise comparison matrix 
indicates how much more important criterion i (say. Topography) is than criterion j 
(say. Operation and Maintenance cost). Thus, ap 2 = 5 indicates that Topography is 
strongly more important than Operation and Maintenance Cost. 



ri 5 9 3 5-1 



5 



1 5 



3 3 




^3511 

43711 
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Fig. 3.1 Hierarchical AHP 
model 




Following the two-step procedure (Step 2), we obtain the weights for each criteria: 
Each entry in column i of A is divided by the sum of the entries in column i, which 
yields the matrix Anorm- Next, W is estimated as the average of the entries in row i 
of Afiorm- 



IF = [ wi W2 W3 W4 ws ] = [0.50 0.10 0.03 0.19 0.19] 

Applying Eq. 3.1 in order to check for consistency, the Eigen value obtained is 
^max — 5.29 and the Consistence Ratio is 0.065, which is less than the allowed value 
of 0.1. Thus, there is a good consistency in the judgments made. 

With respect to Topography, Operation and Maintenance cost, Land use. Infrastruc- 
ture and Investment cost, the following comparison matrices show how well each 
alternative ‘scores’ for each objective: 



Topography Oper. and maint. costs Land use Infrastructure Investment cost 



I 


6 


3 


4 ■ 




'1 


2 


4 


6 " 




'1 


3 


5 


1 n 
6 




'1 


3 


6 


1 n 
5 




'1 


3 


1 

5 


1 -| 

2 


1 

6 


1 


6 


1 

2 




1 

2 


I 


2 


3 




3 


I 


6 


1 

2 




3 


I 


3 


1 

2 




3 


1 


1 

3 


3 


1 

3 


6 


1 


3 




1 

4 


1 

2 


1 


2 




5 


6 


1 


3 




6 


3 


1 


1 

2 




5 


3 


1 


4 


1 

L 4 


2 


.3 


1 _ 




1 


1 


1 

2 


1 _ 




.6 


2 




1 _ 




.5 


2 


2 


1 _ 




_2 


3 


1 

4 


1 _ 



Table 3.4 shows the weight vector for each alternative obtained from the compar- 
ison matrices and Table 3.5 shows the Eigen value (kmax) and the Consistence Ratio 
for each criterion. All the Consistence Ratios are less than the allowed value of 0. 1 . 
Thus, there is a good consistency in the judgments made. 

Eollowing step 4, we now synthesize the (criteria) weights with the score of each 
alternative on each criteria to obtain an overall score for each alternative (Table 3.6). 
Applying Eq. 3.3 La Bragma is the best ranked alternative, followed by Cilda, San 
Pedro, with the Estacas alternative being the least preferred alternative. 
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Table 3.4 Weight vectors 



Cl 


Cl 


C3 


C4 


C5 


■ 0.527" 
0.064 
0.291 
-0.117. 




"0.519" 

0.260 

0.140 

.0.081. 




"0.064" 

0.128 

0.552 

.0.256. 




"0.064" 

0.170 

0.347 

.0.419. 




"0.083" 

0.257 

0.532 

.0.128. 


Table 3.5 Eigen value and consistence ratio for each criteria 




Cl 




Cl 




C3 


C4 


C5 


^max 

Cl 

Rl 


4.12 

0.046 


4.01 

0.004 




4.21 

0.079 


4.14 

0.051 


4.11 

0.040 



Table 3.6 Overall score for 
each alternative 



Alternative 


Score 


A 1 La Bragia'a 


0.343 


Al Cilda 


0.341 


A 3 San Pedro 


0.176 


A 4 Estacas 


0.141 



3.3 Concluding Remarks 

AHP is a decision-making tool that works well when the available data for the 
decision-makers are difficult to quantify or when it is necessary to rely on expert 
opinions or preferences in some intangible aspects. However, the method has been 
frequently criticized for its inability to adequately accommodate the inherent uncer- 
tainty and imprecision associated with certain environments. In addition, it is difficult 
to subjectively scale a concrete quantitatively number to enable the pairwise compar- 
ison without losing some degree of accuracy. Despite these drawbacks, the relative 
ease with which it handles multiple criteria and performs qualitative and quantitative 
data has favored its use as a decision-making method. 

AHP has been widely used in various domains related to the renewable energies. 
Ramanathan and Ganesh [2] performed a multiple objectives analysis to evaluate 
energy alternatives for lighting in households by employing goal programing and 
AHP. Mohsen and Akash [3] evaluated solar water heating systems in comparison 
with other conventional systems in Jordan and Akash et al. [4] performed a com- 
parison between the different electricity power production options in Jordan. Wang 
and Feng [5] discussed the relationship between China’s rural energy and its sustain- 
able development calculating the weighting of an index system for appraising the 
sustainable development of rural energy using AHP. Nigim et al. [6] used AHP to 
assist communities in the pre-feasibility ranking of alternative RE projects to inte- 
grate the use of RE sources with conventional fuel sources in order to minimize 
the dependence on imported energy sources and the negative environmental impact 
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of generating energy. Aras et al. [7] determined the most convenient location for a 
wind farm observation station to be built on the campus of a university. Kablam [8] 
presented an AHP-based modeling framework for the prioritization of energy conser- 
vation policy instruments in Jordan and Vashishtha and Ramachandran [9] evaluated 
demand side management strategies for reducing capacity shortages and improving 
system load factors in the Indian power sector from a multiobjective perspective 
using AHR Nagesha and Balachandra [10] identified and prioritized relevant bar- 
riers to energy efficiency in small-scale industry clusters based on the perceptions 
and experience of entrepreneurs and Kon Lee et al. [11] determined the priorities in 
technology development for the energy efficiency and greenhouse gas control plans 
in Korea. Quintero et al. [12] performed a comparative analysis of the technological 
options using different feedstocks for fuel ethanol production and Kon Lee et al. 
[13] evaluated the potential of Korea to be competitive in development of hydrogen 
energy technology. Lee et al. [14] proposed a MCDM model based on AHP associ- 
ated with benefits, opportunities, costs and risks to help select a suitable wind farm 
project in China. 
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Chapter 4 

Weighted Sum Method and Weighted 
Product Method 



In this chapter we look at two simple multi-criteria decision-making methods, the 
Weighted Sum method and the Weighted Product method. In the Weighted Sum 
method the score of an alternative is equal to the weighted sum of its evaluation 
ratings, where the weights are the importance weights associated with each attribute. 
In the Weighted Product method, instead of calculating sub-scores by multiplying 
performance scores times attribute importance, performance scores are raised to the 
power of the attribute importance weight. 

4.1 Weighted Sum Method 

The Weighted Sum method is a straightforward method, especially used in single 
dimensional problems. If there are m alternatives and n criteria, then the best alter- 
native is the one that satisfies the following expression: 

./ 

Ksm^ Max^aijWj ( 4 . 1 ) 

i 

for i — 1, 2, . . . , m where is the weighted sum method score of the best 

alternative, n is the number of decision criteria, is the actual value of the ith 
alternative in terms of the 7 th criterion and Wj is the weight of importance of the 7 th 
criterion. The total value of each alternative is equal to the sum of products. Difficult 
with this method merges when it is applied to multi-dimensional decision-making 
problems. In combining different dimensions, and consequently different units, the 
additive utility assumption is violated [1]. 

4.2 Weighted Product Method 

The Weighted Product method is similar to the Weighted Sum method. The main 
difference is that instead of addition in this model there is multiplication. Each 
alternative is compared with the others by multiplying a number of ratios, one for 
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4 Weighted Sum Method and Weighted Product Method 



Table 4.1 Evaluation matrix 



Alternative Power {Gw h/year) NPV{Euros * 10°) Maturity(l-5) Emissions(fC 02 />’) Jobs 



Ai Energetic 
cultivations 


1,709 


7.50 


2 


2.25 


15 


A 2 Eorest and 
agricultural 
wastes 


1,815 


8.65 


5 


2.65 


22 


A 3 Eanning 
industrial 
wastes 


698 


9.30 


4 


2.46 


18 


A 4 Eorest 
industrial 
wastes 


698 


9.65 


3 


3.10 


25 



each criterion. Each ratio is raised to the power equivalent to the relative weight of 
the corresponding criterion. In general, in order to compare the alternatives Ak and 
Ai the following product is obtained: 



R 




(4.2) 



where n is the number of criteria, atj is the actual value of the ith alternative in terms 
of the y'th criterion and Wj is the weight of importance of the yth criterion. If R(^) 
is greater than one, then alternative is more desirable than alternative A/ (in the 
maximization case). The best alternative is the one that is better than or at least equal 
to the other alternatives [2]. 



4.3 Application 

In this section, to show the use of the Weighted Sum and Weighted Product methods, 
the selection among four different Biomass plants for electric generation shown in 
Table 4. 1 . is considered. The alternatives will be evaluated according to the follow- 
ing criteria: Electric generation (Guih/year), Net Present Value (NPV), Maturity 
(M), tons of emissions of CO 2 avoided per year (tC02ly) and number of Jobs 
created (I). 

By applying Eq. (4.1) the weighted sum method scores of the four alternatives are 
calculated. Erom Table 4.2. alternative A4, a Biomass plant with forest and industrial 
wastes is considered the best ranked alternative by the Weighted sum method. 

In order to compare the alternatives by the Weighted Product Method, Eq. (4.2) 
is applied. The results, shown in Table 4.3, indicate that alternative A2, a Biomass 
plant with forest and agricultural wastes with R= 1.60, is now the best alternative. 
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Table 4.2 Weighted sum 
method results 



Alternative 


At 


2,475,173 


A2 


2,860,186 


^3 


3,036,073 


A4 


3,205,074 



Table 4.3 Weighted product 
method results 



Alternative 


Al 


Az 


A3 


A4 


Al 


1 


0.63 


1.20 


1.15 


Az 


1.60 


1 


1.21 


1.27 


^3 


1.32 


0.83 


1 


1.05 


A 4 


1.26 


0.79 


0.95 


1 



4.4 Concluding Remarks 

Applications of the Weighted Sum method to sustainable energy systems can be found 
in Jovanovic et al [3], Afgan and Carvalho [4, 5, 6], Liposcak et al. [7], Pilavachi et 
al [8], Begic and Afgan [9], Afgan et al. [10], Patlitzianas et al [11] and Renn [12]. 
Due to the simplicity of the Weighted Sum Method, it often lies behind numerous 
other much more complex aggregation methods such as AHP. With regard to the 
Weighted Product Method a major drawback is that it over-values extremes leading 
to undesirable results because it considerably favors/disfavors the final evaluation of 
any alternative which, relative to one criterion, is far from the average. 
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Chapter 5 

PROMETHEE 



PROMETHEE (Preference Ranking Organization Method for Enrichment Evalu- 
ation) is a multi-criteria decision-making method developed by Brans et al [1]. 
By 1994, the method has been extended to encompass six ranking formats: 
PROMETHEE I (partial ranking), PROMETHEE II (complete ranking), PROME- 
THEE III (ranking based on intervals), PROMETHEE IV (continuous case), PROM- 
ETHEE V (net flows and integer linear programming) and PROMETHEE VI (rep- 
resentation of human brain). The method uses the outranking methodology to rank 
the alternatives combined with the ease of use and decreased complexity. Based on 
extensions of the notion of criterion, the method is well adapted to problem where a 
finite number of alternative actions are to be ranked considering several criteria. 

5.1 PROMETHEE I and II 

The PROMETHEE method is implemented in five steps. In the first step, a preference 
function showing the preference of the decision-maker for an action a with regard to 
another action b, will be defined separately. The second step concerns the comparison 
of the suggested alternatives in pairs to the preference function. As a third step, the 
outcomes of these comparison are presented in an evaluation matrix as the estimated 
values of every criterion for every alternative. The ranking is realized in the two final 
steps: the fourth step includes the PROMETHEE I method application for partial 
ranking and afterward the fifth step includes PROMETHEE II method for complete 
ranking of the alternatives. 

Step 1. Define the Preference function 

Given the preference of the decision-maker for an action a with regard to of a 
set of possible actions K, the preference function, which is a function of the differ- 
ence between two alternatives for any criterion, will be defined separately for each 
criterion; its value will be between 0 and 1. The smaller the function, the greater 
the indifference of the decision-maker; the closer to 1, the greater his preference. 
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5 PROMETHEE 



Fig. 5.1 Type I criterion 



p(x) 

1 ~ 



0 



► 



X 



In case of strict preference, the preference function will be 1. The associated prefer- 
ence function P{a, b) of a with regard to b will be defined as [2]: 



P{a, b) 



0, for f(a) < f(b) 

p[fia), f(b)l for f(a) > f(b) 



(5.1) 



For concrete cases, it seems reasonable to choose p{-) functions of the following 
type 



Plfia), fib)] = p[f(a) - fib)] (5.2) 



depending on the difference between the values/ (a) and/(£>). Six types of functions 
cover most of the cases occurring in practical applications, namely usual criterion, 
quasi-criterion, criterion with lineal preference, level criterion, criterion with linear 
preference and indifference area and Gaussian criterion. In order to indicate clearly 
the areas of indifference in the neighbourhood of/(a) and/(l;), we write: 



.T = fia) - fib) 



and the function Plix) is graphically represented so that: 



//(x) = 



Pia, b), X > 0 
Pib, a), X <0 



1. Type I: Usual Criterion. 
Let pix) he 



pix) = 



0, for Vx < 0 

1, forVx>0 



(5.3) 



(5.4) 



(5.5) 



In this case, there is indifference between a and b only when fia) — fib). 
As soon as these values are different the decision-maker has a strict preference 
for the action having the greatest value. His preference function equals 1 and 
Hix) is given in Fig. 5.1. If the decision-maker identifies the criterion /(.) as 
being of type I, no particular parameter has to be defined. 
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Fig. 5.2 Type II criterion 



Fig. 5.3 Type III criterion 



p( 




i 


k 








1 

1 

1 

1 

1 

1 

1 

^ ► 



-10 1 




2. Type II: Quasi-Criterion. 
In this case. 



p{x) = 



0 , 

1 , 



for X < I 
for X > I 



(5.6) 



and for the particular criterion, a and b are indifferent as long as the difference 
between/(fl) and/(£>) does not exceed Z; if not, the preference becomes strict. 
When the decision-maker identifies the criterion /(•) as being of Type II, only 
the parameter / has to be defined. The function H(x) is given in Fig. 5.2. 

3. Type III: Criterion with Linear Preference 
Let p{x) he 



p(x) = 



— for X < m 

m, — 

1 , for X > m 



(5.7) 



Such an extension of the notion of criterion allows the decision-maker to prefer 
progressively ato b for progressively larger deviations between f{a) and f(b). 
In this case H(x) is given in Fig. 5.3. The intensity of preference increases 
linearly until this deviation equals m, after this value the preference is strict. 
If the decision-maker considers that a particular criterion is of Type III, he has to 
define only the value m from which strict preference is considered. 

4. Type IV: Level-Criterion 



p(x) = 



0 , 

_j_ 

2, 

1 , 



for x < q 

for q < x <q + p 

for x > q + p 



(5.8) 



In this case, a and b are considered as indifferent when the deviation between 
f{a) and f(b) does not exceed q. Between q and q + p the preference is weak 
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Fig. 5.4 Type IV criterion 



Fig. 5.5 Type V criterion 



P(x) 





1 
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(1/2), after this value the preference becomes strict. H{x) has the form depicted 
in Fig. 5.4. The decision-maker can easily fix q and p when he feels that the 
particular criterion /(•) is of Type IV. Criteria with more than two levels can also 
be considered, for instance, when several given norms seem relevant. 

5. Type V: Criterion with Linear preference and Indifference Area. 

This time we consider for pix): 



p{x) = 



0, for X < 

for s < X < -h r 

1 , for X > -h r 



(5.9) 



In this case the decision-maker considers that a and b are completely indifferent 
as long as the deviation between f{a) and f(b) does not exceed s. Above this 
value the preference grows progressively until this deviation equals s + r. H{x) 
is then given by Fig. 5.5. Two parameters, s and r, have to be defined when a 
particular criterion has to be identified as being of this type. 

6. Type VI: Gaussian Criterion 
Let p{x) be: 



pix) = 



0 , 

_ 

1 - e 



for X < 0 
for X > 0 



(5.10) 



If a particular criterion is of the Gaussian type, the preference of the decision- 
maker still grows with the deviation x. The value of a is the distance between 
the origin and the point of inflexion of the curve. H{x) then is given by Fig. 5.6. 
In this particular case only the value of a has to be dehned by the decision-maker. 
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Fig. 5.6 Type VI criterion 



P(x) 




Step 2. Calculate the Preference index 

Suppose every criterion has been identified as being one of the six types considered 
so that the preference functions Pu {a , b) have been defined for each h — 1 , 2 ,..., A:. 
Next, for each couple of actions a, b e K , v/e first define a preference index for a 
with regard to b over all the criteria. Let 

1 ^ 

7T{a,b) = -^Ph(a,b) (5.11) 

* h=i 

be such a preference index, which gives a measure of the preference of a over b for 
all the criteria. The closer to 1, the greater the preference. 

Step 3. Construct the valued outranking graph 

The values obtained in Step 2 determine the valued outranking graph, the nodes of 
which are the actions of K, so that for all a, A; e K, the arc (a,b) has the value 
jt(a,b). Let us define, for each node in this valued outranking graph, the outgoing 
flow: 

(/)+(«) = ^ :r(fl, x) (5.12) 

xeK 

and the incoming flow: 

<p~ {a) (5.13) 

x(iK 



The larger (p~^(a), the more a dominates the other actions of K. The smaller 
4>~(a), the less a is dominated. 

Step 4. Ranking the actions hy a partial preorder 

If the decision-maker wants to rank the actions of K from the best to the weakest 
one, the problem consists in using the outranking graph to build a total preorder on 
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K, or a partial one. Let us define the two total preorders and (P , I ) 

such that 



a P"*" b 


if (p'^{a) > (p'^ib) 


(5.14) 


a P~ b 


if (p~{a) < (p~{b) 


a b 


if (p^(a) — (p'^ib) 


(5.15) 


a I~ b 


if (p~{a) — <p~{b) 



We then obtain the following partial preorder R) by considering their 

intersection: 



a outranks b (a P^^^ b) 

a is indifferent to b (a b) 
a and b are incomparable {a R b) 



a b and a P~ b 
if •! a P"*" b and a I~ b 
a b and a P"*" b 
if a b and a I~ b 
otherwise 



(5.16) 



This is the PROMETHEE I partial relation. It offers the decision-maker a graph in 
which some actions are comparable, while some others are not. 



Step 5. Ranking the actions by a total preorder 

Suppose a total preorder (complete ranking without incomparabilities) has been 
requested by a decision-maker. We then consider for each criterion a e K the net 
flow: 



(/)(a) = (p~^{a) — (j) (a) (5.17) 

which can easily be used for ranking the actions: 

a outranks b {a P^^^ b) if <p(a) > 4>(b) 

a is indifferent to b {a 1^^^ b) if 4>{a) — <p{b) 

This is the PROMETHEE II complete relation. All the actions of K are now 
completely ranked but this relation is poorer in information and less realistic because 
of the balancing effects between outgoing and incoming flows. 



5.2 Application 

Let us consider the multi-criteria problem proposed by Goumas and Lygerou [3] con- 
sisting in the selection of four alternative exploitation scenarios for the low enthalpy 
geothermal field in a rural community located in Northern Greece. Each exploita- 
tion scenario includes a number of possible uses such as: greenhouse heating, sub- 
soil heating, drying of agricultural products and water heating for fish farming. For 
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Table 5.1 Alternatives for a geothermal field [3] 



Alternatives 


NPV (lO^GRD) 


Jobs (N) 


Energy used (10*^ cal/year) 


Risk 




Type II 


Type III 


Type V 


Type IV 




1=100 


m=I0 


s= 1 


q=I 








r=3 


p = 3 


Ai 


373 


44 


22.4 


3 


A2 


3706 


99 


51.0 


6 


A3 


3809 


100 


49.3 


8 


A4 


3860 


104 


47.5 


9 


weight 


0.30 


0.30 


O.IO 


0.30 



more details on the formation of alternative exploitation schemes see Goumas and 
Lygerou [3], four criteria have heen selected for the evaluation: the Net Present Value 
of the investment (NPV in Greek Drachmes), the creation of new jobs (Jobs), the 
Energy consumed which indicates the degree of exploitation of the resource and the 
Risk associated either with new products and/or the application of new technology 
or with uncertainties for future changes in the economic parameters. The index is 
expressed in a scale of 1 (no risk) to 20 (very risky). The criteria NPV, and Jobs have 
to be maximized whereas Energy used and Risk have to be maximized. 

The data for the evaluation, the type of each criterion and the the values of the 
corresponding parameters are shown in Table 5.1. In Goumas and Lygerou [3], 
the generalized criterion with linear preference was applied in all criteria. To show 
the use of different criteria, the following types of criterion have been considered: 



1. Net Present Value. Type of criterion II: Quasi-criterion. 



p(x) = 



0, for X < 100 

1, for X > 100 



In this case two alternatives are indifferent as long as the difference between them 
regarding the investment required does not exceed 1 00 GRD . If not, the preference 
becomes strict. 

2. Creation of new Jobs. Type of criterion III. Criterion with Linear Preference 



p{x) = 



for X < 10 
1, for X > 10 



This criterion allows the decision-maker to prefer progressively alternative a to 
b for progressively larger deviations between f{a) and f(b). The intensity of 
preference increases linearly until this deviation equals 10 jobs. After this value 
the preference is strict. 

3. Energy consumed. Type of criterion V: Criterion with Linear Preference and 
Indifference Area 
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Table 5.2 Values of 
n(Ai,Ak) 



Table 5.3 Data for preorders 
and 





Ai 


^2 


^3 


A 4 


Ai 


- 


0.375 


0.500 


0.500 


^2 


0.500 


- 


0.125 


0.125 


^3 


0.500 


0.332 


- 


0 


A 4 


0.500 


0.707 


0.315 


- 




Ai 


^2 


^3 


A 4 




1.375 


0.750 


0.832 


1.522 


<p-{Ak) 


1.500 


1.414 


0.940 


0.625 



p{x) = 



0 , 

x-l 

3, 

1 , 



for X < 1 
for 1 < X < 4 
for X > 4 



In this case the decision-maker considers that two alternatives are completely 
indifferent as long as the deviation between them does not exceed 1 * lO'^ cal/year. 
Above this value the preference grows progressively until this deviation equals 
4 * 10^^ cal/year hours. 

4. Risk. Type of criterion IV: Level Criterion 



p(x) = 



0, for X < 1 

Y for 1 < X < 4 

1 , for X > 4 



In this case, two alternatives are considered as indifferent when the deviation 
between them does not exceed 1. Between 1 and 4 the preference is weak, and 
after this value the preference becomes strict. 

Using Eqs. 5.1, 5.3 and 5.1 1 in the right way according to whether the criterion 
has to be minimized or maximized, we obtain the preference index for each couple 
of alternatives. Thus, for alternatives (Ai) and (A 2 ) the preference index will be: 

0(^1, M) — M) + PiiAi, Aj) + P 3 (A[, A 2 ) + Pa{A\, A 2 )] 

= ^[0 + 0+ 1 + = 0.375 

Table 5.2 shows the values 7 t(A, , A/;) for every couple of actions. These values 
determine the outranking graph, each arc (A;, Ak) having the value 4>(Ai, At). 

Let us first suppose that a partial relation would be useful to the decision-maker. 
We therefore apply the PROMETHEE I technique. According to Eqs. 5.12 and 5.13 
we complete Table 5.3 with the outgoing and the incoming flow. 

Nexf, we obfain the preorders and P~ from Table 5.3, the intersection of 
which is: 



A4P^‘)Ai, AaP^^^Aj, AaP^^^A^, A^P^^'^Aj 
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Table5.4 Net flows 



A I A2 A3 A4 

0(A,) -0.125 -0.644 -0.108 0.897 



Supposing now that the decision-maker requests a total preorder, we can then 
use PROMETHEE II. According to Eq. 5.17 we calculate the net flows shown in 
Table 5.4. Erom Table 5.4, alternative A 4 is considered the best ranked alternative. 



5.3 Concluding Remarks 

PROMETHEE methods are well adapted to problems where a finite number of alter- 
native actions are to be ranked considering several criteria. The methods are rather 
popular in the world of outranking methods and one of the reasons for this popular- 
ity is the existence of user friendly softwares such as PROMCALC. However, when 
applying the method several considerations must be taken into account. For example, 
the decision-maker can express his preferences between two actions on all criteria 
and the importance he attaches to the criteria on a ration scale, for all criteria the 
difference between evaluations must be meaningful. 

Mladineo et al [4] used PROMETHEE method to deal with the problem of small 
hydropower construction in a wide area covering many regions in order to make the 
units costs as low as possible. Georgopoulou et al. [5] designed and implemented 
a decision support system by means of PROMETHEE II for sustainable renewable 
energy exploitation in Greece. Diakoulaki and Karangelis [ 6 ], Goumas et al. [7] 
and Goumas and Lygerou [3] applied the method to geothermal energy projects 
and Haralambopoulos and Polatidis [ 8 ] to the exploitation of a geothermal resource 
in the island of Chios (Greece). Pohekar and Ramachandran [9] used the method 
for promoting parabolic solar cooker in India. Madlener et al. [10] explored future 
renewable energy options in Austria and Tsoutos et al. [11] the sustainable energy 
planning on the island of Crete. 
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Chapter 6 

ELECTRE 



The methods belonging to the ELECTRE family (ELimination Et Choix Traduisant 
la REalite) devised by Roy [1] grew from the idea that rigorous mathematical axioms 
are unsuited to describe a complex reality such as the one of a contradiction-laden 
decision process. ELECTRE methods are based on the partial aggregation of prefer- 
ences, which is based on the principle that the decision-maker is not perfectly rational 
and can therefore express, for each pair of actions (a and b), not only the preference 
(P) or indifference (I) toward one action or the other, but also a less marked preference 
labeled ‘weak’ (Q) or a situation of incomparability (R). 

ELECTRE methods are not compensatory, a high degree of dissatisfaction regard- 
ing the performance of one alternative compared to another for one specific criterion 
may not necessarily be compensated for by good performance with regard to the other 
remaining criteria. It may therefore happen that the alternatives differ so widely from 
one another that the decision-maker is unable to express a preference as they are not 
comparable. 



6.1 ELECTRE III 

The distinctive characteristic of ELECTRE III is its ability to handle a data-set 
affected by a high degree of uncertainty. It also has the advantage of being less 
sensitive to any changes in data, which explains why ELECTRE III is very well-suited 
to environmental applications and it also involves some aspects neglected by other 
methods [2]. The method is capable of handling discrete criteria both quantitative 
and qualitative in nature. 

The evaluation procedures of the ELECTRE III method encompass the estab- 
lishment of the threshold function, disclosure of concord index and discord index, 
outranking degree and the ranking of alternatives [3]. 
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Step 1. Establishment of the threshold function 

Let qj and pj represent the indifference threshold and preference threshold, respec- 
tively, with regard to the y'th criterion. Indifference and preference 
thresholds characterize the acceptance of imprecision in the judgment by consid- 
ering two entities as indifferent when their individual performance in each criterion 
j differs less than a specified amount qj. In particular, qj indicates the minimum 
boundary of uncertainty associated with the performed calculations, while p j can 
be considered as the maximum boundary of error, connected to the performed cal- 
culations. Moreover, the transition from indifference to preference is made gradual, 
changing linearly from qj to pj. Threshold functions are linked to the idea of a 
pseudo-criterion. A criterion gj is a pseudo-criterion if there exists two threshold 
functions q{g) and p{g) that represent the indifference and preference threshold, 
respectively, so that for g{a) > gib) there are situations in-between in which a weak 
preference for option a over b denoted aQb applies as follows [3]. 

1. gia) > gib) + pigib)) ^ aPb 

2. gib) +qigib)) < gia) < gib) + pigib)) aQb 
3- gib) < gia) < gib) + qigib)) alb 

where P refers to strong preference, Q weak preference, I indifference and gia) and 
gib) are the criteria value of alternative a and b, respectively. 

The establishment of a threshold function has to satisfy the subsequent constraint 
equations: 

1. gia) > gib) gia) +qigia)) > gib) + qigib)) 



where gjia) is the evaluation value of alternative a on criterion j. The threshold 
values may be constant (a is determined and /I = 0) or proportional ifi > 0). 

Step 2. Concordance index and Discordance index 

A Concordance index C ia,b) is computed for each pair of alternatives: 



gia) + pigia)) > gib) + pigib)) 



2. for all criteria, pig) > qig) 

Pjigjia)) and qjigjia)) can be assessed according to Roy’s formula [4]: 



Pjigjia)) = Q!p + Ppgjia) 



qjigjia)) = + Pqgjia) 




( 6 . 1 ) 
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where kj is the weight of criterion j and cj (a , b) is the outranking degree of alternative 
a and alternative b under criterion j, calculated as follows: 



Cj{a, b) 



1 

0 

Pj+9j(a)-9j(b) 

Pj-9j 



if gj {a) + qj > gj(b) 
if gj{a) + Pj < gj{b) 
otherwise 



A discordance index dj (a, b) for each criterion is then defined as: 



dj{a, b) 



0 

1 

9j(b)-gj{a)-pj 

'’J-PJ 



ifgjia) + Pj > gj{b) 
if gj {a) + Vj < gj(b) 
otherwise 



( 6 . 2 ) 



(6.3) 



where Vj is the veto threshold value defined for each criterion j. The veto threshold 
value is the limit value of the difference gj (b) — gj (a), over which it is reasonable to 
reject the hypothesis of outranking of a over b, with regard to the considered criterion. 
The veto threshold is aimed at capturing situations in which very bad scores in any 
criterion should prevent an entity of being classified in the best category, or if these 
bad scores should force it to be classified in the worst category independent of having 
very good scores in all other criteria. 



Step 3. Outranking degree 

The next step is to combine the concordance and discordance measure for each 
pairs of alternatives to produce a measure of the degree of outranking; that is, 
a credibility matrix which assesses the strength of the assertion a is at least as good 
as b. The credibility degree for each pair of alternatives is defined by S{a, by. 



S{a,b) 



C{a,b), if dj{a,b) < C(a,byi j 

C(fl,&)*n/e 7 (fl,fo) otherwise (6.4) 

^(i - dj{a,b))l'{i - cj(a,b)). 



where J(a,b) is the set of criteria for which dj{a, b) > C{a, b) 

This formula assumes that if the strength of the concordance exceeds that of the 
discordance, then the concordance value should be modified. Otherwise, we are force 
to question the assertion that aSb and modify C(a, b) according to the above equation. 



Step 4. Ranking of alternatives 

To obtain a complete ranking of the alternatives, the normal ranking method of ELEC- 
TRE III uses a structured algorithm via two intermediate ranking procedures: one 
descending where alternatives are classified from the best to the worst (descending 
distillation) and the other ascending from the worst to the best (ascending distillation) . 
The partial rankings obtained take into account indifference but not incomparabil- 
ity and they are therefore intertwined in order to obtain the final ranking in which 
incomparability is also taken into account. Eollowing Li and Wang [5] the weakness 
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of the normal ranking of ELECTRE III is that it requires an additional threshold to be 
introduced and the ranking of the alternatives depends on the size of this threshold for 
which there exists no ‘correct’ value. Aiming at the ranking problems in ELECTRE 
III, Li and Wang [5] present a new ranking method by introducing three concepts, 
the concordance credibility degree, the discordance credibility degree and the net 
credibility degree. 

1 . The concordance credibility degree is defined by 

d)+(x/) = ^ Sixi,xj), Vxi e Z (6.5) 



The concordance credibility degree is the measure of the outranking character of 
Xi (how X; dominates all the other alternatives of X). 

2. The discordance credibility degree is defined by 

<t>“(x,) = S(xj,xi), Vx/ e X (6.6) 

Xj£X 



The discordance credibility degree gives the outranked character of x,- (how x; is 
dominated by all the other alternatives of X). 

3. The net credibility degree is defined by 

0 (X; ) = 0+ (xi ) - 0- (Xi ) , Vxi e Z (6.7) 

The net credibility degree represents a value function, where a higher value reflects 
a higher attractiveness of alternative x/ . 

Next, all the alternatives can be completely ranked by the net credibility degree. 



6.2 ELECTRE IV 



The basic difference between ELECTRE III and IV is that no weights for the cri- 
teria are introduced in ELECTRE IV, which does not mean that the weights of the 
criteria are considered to be equal. The ELECTRE IV method uses the following 
notation [3]; 

• mp(a, b): the number of criteria for which alternative a is strictly preferred to 
alternative b. 

• niqia, b): the number of criteria for which alternative a is weakly preferred to 
alternative b. 

• mi{a,b)'. the number of criteria for which alternative a is considered indifferent 
to alternative b but such that a has a better criterion value than b. 

• nioia, b) — Moib, a): the number of equal criterion values of alternative a and 
alternative b. 
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Table 6.1 Evaluation matrix 



Alternatives 


Power (Mw) 


NPV(Euros * 10*’) 


Maturity (1-5) 


Emissions (tC02ly) 


Jobs 


Ai Energetic 
cultiv. 


245 


7,500,000 


2 


2,250,000 


15 


A 2 Eorest and 
agric. wastes 


260 


8,650,000 


5 


2,650,000 


22 


A 3 Eanning 
indust, wastes 


100 


9,300,000 


4 


2,460,000 


18 


A 4 Eorest indust, 
wastes 


230 


9,650,000 


3 


3,100,000 


25 



The outranking relations in the ELECTRE IV method are the following: 

1 . Quasi-dominance Sq 

aSqb nip{b, a) + niq{a, b) — 0 and rrii(b, a) < 1 -t- nii{a, b) -\- niq{a, b) + 
mp{a, b) 

2. Canonic-dominance Sc 

aScb 4^ mp(b, a) — 0 and mq(b, a) < mp{a, b) and mq{b, a) + m; {b, a) < 
1 4- nii{a, b) + ntq{a, b) + nip(a, b) 

3. Pseudo-dominance Sp 

aSpb 44 mp(b, a) — 0 and mq{b, a) < mq(a, b) + nip{a, b) 

4. Veto-dominance 

aSyb 44 rripib, a) = 0 if mp(b, a) = 1 and no bPVja, V j and mp(a, b) > '-j in 
which bPVja 44 gjib) > gjia) + Vj(gjia)) 

The partial preorder is constructed in a similar way as in ELECTRE III. 



6.3 Application 

The selection among four different Biomass plants for electric generation is con- 
sidered in order to show the use of ELECTRE III method. The evaluation matrix 
is shown in Table 6.1. The alternatives will be evaluated according to the following 
criteria: Power (Mw), Net Present Value (NPV), Maturity (M), tons of emissions of 
CO 2 avoided per year (tCOjIy) and number of Jobs created (7). 

Once the criteria has been defined, it is necessary to assign weights to these crite- 
ria. In ELECTRE methods, weights are perceived as true coefficients of importance 
assigned to the criteria, and not just as technical devices for translating the perfor- 
mance in the criteria considered into a common value measure. Therefore, they are 
scale-independent, that is, they are not linked to the scales in which each criterion 
is measured, thus making them easier to be specified by managers [6]. The most 
used technique when determining the weights in ELECTRE methods is the Sismos 
approach [7]. Table 6.2 shows the weighting procedure, according to this technique. 
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Table 6.3 Weights and threshold values 



Weights 


Power 

0.23 


NPV 

0.23 


Maturity 

0.10 


Emissions 

0.34 


Jobs 

0.10 


Indifference threshold (q) 


15 


300,000 


1 


250,000 


3 


Preference threshold (p) 


30 


500,000 


2 


400,000 


5 


Veto threshold (u ) 


70 


750,000 


3 


500,000 





based on a first ranking of criteria and a subsequent assignment of weights depending 
on each criterion rank. 

Once the problem has been defined, the alternatives generated and the weights 
assigned, threshold values must be considered. Procedures for choosing appropri- 
ate threshold values are addressed by Rogers and Bruen [ 8 ]. In our case, we have 
considered the indifference, preference and veto threshold shown in Table 6.3. 

Next, the concordance index is computed. The concordance calculations for alter- 
natives Ai and A 2 are as follows: 



Cl (A 3 , A 4 ) 
C 2 (A 3 , A 4 ) 

C 3 (A 3 , A 4 ) 
C 4 (A 3 , A 4 ) 

C 5 (A 3 , A 4 ) 



0 since 100 -|- 15 < 230 
and 100 -h 30 < 230 

0.75 since 9, 300, 000 -h 300, 000 < 9, 650, 000 
and 9, 300, 000 -h 500, 000 > 9, 650, 000 



then 



500, 000 -h 9, 300, 000 - 9, 650, 000 



500, 000 - 300, 000 
since 4 -|- 1 > 3 

since 42, 460, 000 -h 250, 000 < 3, 100, 000 
and 2, 460, 000 -h 400, 000 < 3, 100, 000 
since 1 8 -|- 3 < 25 
and 18 -I- 3 <25 



= 0.75 



Therefore CfAs, A 4 ) = 0.23 * 0.75 -b 0.10 * 1 = 0.27. 

The value of 0.27 measures the strength of the assertion that A 3 is at least as good 
as A 4 . Table 6.4 presents the complete concordance matrix. 

In the concordance matrix we have a measure of the extent to which we are in 
harmony with the assertion that alternative a is at least as good as alternative b. Next, 
to calculate if there is any discordance associated with the assertion aSb the veto 
threshold is used. This veto threshold allows for the possibility of aSb to be refused 
totally if, for any one criterion j, gjia) + vj < gjib). 

Consider criterion C 2 , with a veto threshold of 750,000. Comparing Ai and A 4 , 
it is clear that 7, 500, 000 + 500, 000 < 800, OOOor gj{a) + vj < gj{b). 

Therefore, the discordance index d\{A\, A 4 ) = 1 . Unlike concordance, no aggre- 
gation over criteria takes place; one discordant criterion is sufficient to discard out- 
ranking. 

There now exists a concordance and a discordance measure for each pair of 
alternatives. The next step is to combine these two measures to produce a measure 
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Table 6.4 Concordance 
matrix 



Table 6.5 Credibility matrix 





Al 


A 2 


A 3 


A4 


Al 


1 


0.23 


0.67 


0.33 


A 2 


1 


1 


0.77 


0.43 


^3 


0.77 


0.67 


1 


0.27 


A 4 


1 


0.67 


1 


1 






Al 


A 2 


A 3 


A 4 


Al 


1 


0 


0 


0 


A 2 


1 


1 


0 


0 


^3 


0 


0 


1 


1 


A4 


1 


0.67 


1 


1 



Table 6.6 Concordance, 
discordance and net 
credibility matrix 



Alternative 


4 ,+ 


4> 


0 


Al 


1 


3 


-2 


A 2 


2 


1.67 


0.33 


^3 


1 


2 


-1 


A 4 


3.67 


1 


2.67 



of the degree of outranking, that is, a credibility matrix which assesses the strength 
of the assertion that ‘a is at least as good as V . If the discordance is 1 .00 for any 
criterion, then we have no confidence that aSb and S{a, b) = 0.00. The credibility 
matrix for our multi-criteria problem is shown in Table 6.5. 

The next step in the model is to produce a ranking of alternatives from the cred- 
ibility matrix. Following Li and Wang [5] the concordance, discordance and net 
credibility degree are calculated by using Eqs. 6.5, 6.6 and 6.7 respectively and 
shown in Table 6 . 6 . From Table 6 . 6 , we may conclude that alternative A 4 , Forest 
industrial wastes, is the best ranked alternative. 



6.4 Concluding Remarks 

Applications of ELECTRE methods can be seen in Barda et al. [9] for the location 
of thermal power plants, Georgopoulou et al. [10] for the renewable energy planning 
in a Greek island. Beccali et al. [11] for the diffusion of RE technologies in the 
island of Sardina (Italy), Papadopoulos and Karagiannidis [12] for the optimization 
of decentralized energy systems, Madlener et al. [ 6 ] to assess the performance of 
agricultural biogass plants and Caravallo [2] to assess a selection of production 
processes of thin-film solar technology. 
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Chapter 7 

TOPSIS 



The principle behind the Technique for Order Preference by Similarity to Ideal Solu- 
tions (TOPSIS) is simple. Ideal and a negative-ideal solutions are formed. The ideal 
solution is formed as a composite of the best performance value exhibited by any 
alternative for each attribute and the negative-ideal solution is the composite of the 
worst performance values. The chosen alternative should be as close to the ideal 
solution as possible and as far from the negative-ideal solution as possible. In this 
chapter, the method is applied to the selection of a renewable project for electric 
generation. 



The Technique for Order Preference by Similarity to Ideal Solutions (TOPSIS) 
method is developed by Hwang and Yoon [1] as an alternative to ELECTRE. The 
basic principle is that the selected alternative should have the shortest distance from 
the ideal solution and the farthest distance from the negative-ideal solution in geomet- 
rical sense. The method assumes that each attribute has a monotonically increasing 
or decreasing utility. This makes it easy to locate the ideal and negative-ideal solu- 
tions. Thus, the preference order of alternatives is yielded through comparing the 
Euclidean distances. The TOPSIS method consists of the following steps; 

Step 1. Calculate the normalized decision matrix 

A decision matrix of the m alternatives and n criteria is formulated first. Then the 
normalized value r,/ is calculated as 



where fij is the value of the i th criterion function for the alternative Aj {j — I, . . . , m; 
i — I, ... ,n). 



7.1 TOPSIS 




(7.1) 
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Step 2. Calculate the weighted normalized decision matrix 

The weighted normalized value Vij is calculated as: 

Vij = Winj (7.2) 

where w, is the weight of the i criterion or attribute and X/'=i — 1- 

Step 3. Determine the ideal and negative-ideal solutions 

For benefit criteria the decision-maker wants to have maximum values among the 
alternatives and for cost criteria he wants minimum values among alternatives. Then, 
the ideal solutions (A*) and the negative-ideal solutions (A~) will be: 

A* = . . . , u„*} = limaxjVij \i e I'), (mmjVijli e /")} (7.3) 

A“ = {uj", . . . , u“} = {(min;i;,7|i e I'), ima.XjVij\i e /")} (7.4) 



where I' is associated with benefit criteria, and I” is associated with cost criteria 

Step 4. Calculate the separation measures 

Using the n-dimensional Euclidean distance, the separation of each alternative from 
the ideal solution is given as: 

n 

/ = ! 

Similarly, the separation from the negative-ideal solution is given as: 



^.7 = 



'^{Vij-Vi )2 



(7.6) 



\ ;=i 



Step 5. Calculate the relative closeness to the ideal solution 

The relative closeness of the alternative aj with respect to A* is defined as: 










(D* + DJ) 



(7.7) 



Step 6. Rank the preference order 

Rank the alternatives, sorting by the value C* in decreasing order. Propose as a 
solution the alternative which is the best ranked (maximum) by the measure. 



7.2 Application 
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Table 7.1 Alternatives for electric generation 



Alternatives 


Ai 


Wind power P < 5Mvo 




Wind power 5 < P < lOMni 


A 3 


Wind power 10 < P < 50Mw 


A 4 


Hydroelectric P < lOMiu 


As 


Hydroelectric 10 < P < 25Mw 


^6 


Hydroelectric 25 < P < 5QMw 


Ai 


Solar Thermo-electric P > lOMiu 


As 


Biomass (energetic cultivations) P < 5Mw 


A 9 


Biomass (forest and agricultural wastes) P < 5Mw 


Aio 


Biomass (farming industrial wastes) P < 5Mw 


All 


Biomass (forest industrial wastes) P < 5Mw 


Ai 2 


Biomass (Co-combustion in conventional central) P > 50Mw 


Ai3 


Bio fuels P < 2Mw 



Table 7.2 Criteria used to evaluate the alternatives 





Criteria 


Unit 


Weight 


fl 


Power (P) 


Kw 


0.32 


f2 


Investment ratio (IR) 


Euros/ Kw 


0.09 


h 


Implement period (IP) 


Years 


0.03 


/4 


Operating hours (OH) 


Hours/year 


0.12 


/5 


Useful life (UL) 


Years 


0.13 


/6 


Operation and maintenance costs (OM) 


Euros* 10^^ / Kwh 


0.04 


/7 


Tons of CO 2 avoided (tC02ly) 


Tons 


0.27 



7.2 Application 

Let us consider in this section as a multi-criteria problem, the selection among 13 
different RE alternatives for electric generation shown in Table 7.1. The designed 
systems will be evaluated according to the criteria shown in Table 7.2. The attributes 
considered are: Power (P), Investment Ratio (IR), Implementation Period (IP), Oper- 
ating Hours (OH), Useful Life (UL), Operation and Maintenance costs (OM) and tons 
of emissions of CO 2 avoided per year (tC02jy)- The weights have been estimated 
by AHP The evaluation matrix showing the data considering the 1 3 alternative RE 
projects and 7 selection criteria are shown in Table 7.3. 

Eirstly, the normalized decision matrix is calculated using Eq. 7.1. Then, using 
Eq. 7.2, the construction of the weighted decision matrix, shown in Table 7.4, 
is carried out. Table 7.5 shows the ideal and negative-ideal solutions, determined 
from the weighted normalized decision matrix, taking into account that, for benefit 
criteria the decision-maker wants to have maximum values among the alternatives 
and for cost criteria, he wants minimum values among the alternatives. This is fol- 
lowed by the separation measures of each alternative from the ideal solution and from 
the negative-ideal solution calculated using Eqs. 7.5 and 7.6, respectively. Next, using 
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Table 7.3 Evaluation matrix 



Criteria 


P 


IR 


IP 


OH 


UL 


OM 


tCOily 


Ai 


5,000 


937 


1 


2,350 


20 


1.47 


1,929,936 


Az 


10,000 


937 


1 


2,350 


20 


1.47 


3,216,560 


^3 


25,000 


937 


1 


2,350 


20 


1.51 


9,649,680 




5,000 


1,500 


1.5 


3,100 


25 


1.45 


472,812 


As 


20,000 


700 


2 


2,000 


25 


0.70 


255,490 




35,000 


601 


2.5 


2,000 


25 


0.60 


255,490 


Ay 


50,000 


5,000 


2 


2,596 


25 


4.20 


482,856 


^8 


5,000 


1,803 


1 


7,500 


15 


7.11 


2,524,643 


Ag 


5,000 


1.803 


1 


7,500 


15 


5.42 


2,524,643 


^10 


5,000 


1,803 


1 


7,500 


15 


5.42 


2,524,643 


All 


5,000 


1,803 


1 


7,500 


15 


2.81 


2,524,643 


Ai2 


56,000 


856 


1 


7,500 


20 


4.56 


22,524,643 


Ai3 


2,000 


1,503 


1 


7,000 


20 


2.51 


5,905,270 



Table 7.4 Weighted normalized decision matrix 


Criteria 


P 


IR 


IP 


OH 


UL 


0 and M 


tCOz/v 


Ai 


0.018 


0.012 


0.005 


0.014 


0.037 


0.005 


0.038 


A 2 


0.035 


0.012 


0.005 


0.014 


0.037 


0.005 


0.063 


A3 


0.088 


0.012 


0.005 


0.014 


0.037 


0.005 


0.190 


A 4 


0.018 


0.020 


0.008 


0.019 


0.046 


0.005 


0.009 


As 


0.071 


0.009 


0.010 


0.012 


0.046 


0.002 


0.005 




0.124 


0.008 


0.013 


0.012 


0.046 


0.002 


0.005 


Ai 


0.177 


0.066 


0.010 


0.016 


0.046 


0.014 


0.010 


Ag 


0.018 


0.024 


0.005 


0.045 


0.027 


0.023 


0.050 


Ag 


0.018 


0.024 


0.005 


0.045 


0.027 


0.017 


0.050 


Aio 


0.018 


0.024 


0.005 


0.045 


0.027 


0.017 


0.050 


All 


0.018 


0.024 


0.005 


0.045 


0.027 


0.009 


0.050 


Ai2 


0.198 


0.011 


0.005 


0.045 


0.037 


0.015 


0.095 


Ai3 


0.007 


0.020 


0.008 


0.042 


0.037 


0.008 


0.116 



Eq. 7.7, the relative closeness of every alternative to the ideal solution is calculated. 
The separation measures and the relative closeness to the ideal solution are shown 
in Table 7.6. 

Finally, the ranking of the alternatives sorting hy the value C* is performed. 
The best alternative is one which has the shortest distance to the ideal solution 
and longest distance to the ideal solution. The results shown in Table 7.7 indicate 
that the Biomass plant alternative (Co-combustion in a conventional power plant, 
P > 50Mw) is the best choice, followed by the Windpower (10 < P < 50Mw) 
and Solar Thermoelectric (P > lOMin) alternatives. 



7.3 Concluding Remarks 
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Table 7.5 Ideal (u*) and negative-ideal solutions (v^ ) 



Criteria 





P 


IR 


IP 


OH 


UL 


O and M 


tC02ly 




Max 


Min 


Min 


Max 


Max 


Min 


Max 




0.198 


0.008 


0.005 


0.045 


0.046 


0.002 


0.190 




0.018 


0.066 


0.013 


0.012 


0.027 


0.023 


0.005 



Table 7.6 Separation measures and relative closeness to the ideal solution 





Ai 


^2 


A3 


A4 


A5 


^6 


Ai 


As 


A9 


A 10 


All 


Ai2 


Ai3 


D* 


0.249 


0.221 


0.135 


0.271 


0.248 


0.232 


0.218 


0.242 


0.241 


0.241 


0.241 


0.120 


0.223 




0.066 


0.084 


0.206 


0.054 


0.082 


0.124 


0.161 


0.070 


0.070 


0.070 


0.071 


0.212 


0.126 


q 


0.211 


0.276 


0.605 


0.165 


0.249 


0.349 


0.424 


0.225 


0.226 


0.226 


0.229 


0.639 


0.361 



Table 7.7 Ranking 



Ai2 


A3 


A7 


Ai3 


A6 


A2 


As 


All 


A9 


Aio 


As 


Ai 


A4 


C* 0.639 


0.605 


0.424 


0.361 


0.349 


0.276 


0.249 


0.229 


0.226 


0.226 


0.225 


0.211 


0.165 



7.3 Concluding Remarks 

The main drawback of the TOPSIS method is that it introduces the ranking index 
including the distances from the ideal point and from the negative-ideal point without 
considering their relative importance, an issue that will be solved by using the VIKOR 
method in the following chapter. Applications of the TOPSIS method in the literature 
are very limited. To quote a few examples, the method has been applied to solid waste 
management [2], fuel buses for public transportation [3], sustainability of renewable 
energy options [4], the assessment of offshore wind turbine support structures [5] 
and in uncertain environments, the method has been applied to solar systems and 
energy planning [6, 7]. 
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Chapter 8 

VIKOR 



The Compromise Ranking Method, also known as the VIKOR method (VlseKri- 
terijumska Optimizacija I Kompromisno Resenje in Serbian, means Multicriteria 
Optimization and Compromise Solution), is an effective tool in multi-criteria deci- 
sion making. This method introduces the multi-criteria ranking index based on the 
particular measure of ‘closeness’ to the ‘ideal’ solution, F* . The compromise solu- 
tion Fc is a feasible solution that is the ‘closest’ to the ideal solution, and compromise 
means an agreement established by mutual concessions. In this chapter, the method 
is applied to the selection among 13 different renewable energy options for electric 
generation. 



8.1 VIKOR 

Assuming that each alternative is evaluated according to each criterion function, the 
compromise ranking could be performed by comparing the measure of closeness to 
the ideal alternative and a compromise means an agreement established by mutual 
concessions. The multi-criteria measure for compromise ranking is developed from 
the L/?-metric used as an aggregating function in a compromise programming method 
[1,2]: 



L 



pj 



fij) 

h ift-fr) 




( 8 . 1 ) 



where Lij and Too.y are used to formulate ranking measure; fij is the value of the 
ith criterion function for the alternative Aj\ and n is the number of criteria. The 
compromise ranking algorithm VIKOR is implemented carrying the following steps 

[ 3 ]: 
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Step 1. Determine the best and the worst values of all criterion functions 

The best (/*) and the worst (/”) values of all criterion functions are determined, 
taking into account that, if the ith function represents a benefit then: 



f* = maxjfij 
fr = minjfij 



while if the function represents a cost then: 



( 8 . 2 ) 



A*_ = minjfij 

fi ^maxjfij 



( 8 . 3 ) 



Step 2. Compute the values Sj and Rj 

The values Sj and Rj are calculated by the relations: 






- fij) 

if* - fn 



and 



Rj — maxi 



mif* - fij) 

(ft - fn 



( 8 . 4 ) 



( 8 . 5 ) 



where wt are the weights of criteria, expressing the decision-maker’s preference as 
the relative importance of the criteria. 



Step 3. Compute the values Q j 

The values Q j are calculated by the relation 



Qj = 



S* 



+ (l-t^) 



(Rj - R*) 

(R- - R*) 



( 8 . 6 ) 



maxjRj and v is 



where S* = minjSj', S~ = maxjSj', R* = minjRj\ R~ = 
introduced as a weight for the strategy of maximum group utility, whereas (1 — i>) 
is the weight of the individual regret. The solution obtained by minjSj is with a 
maximum group utility (‘majority’ rule), and the solution obtained by minjRj is 
with a minimum individual regret of the ‘opponent’. 



Step 4. Rank the alternatives 



Rank the alternatives, sorting by the values S, R and Q in decreasing order. The results 
are three ranking lists. 



Step 5. Propose the compromise solution 

Propose as a compromise solution the alternative which is the best ranked by 
the measure Q (minimum), if the following two conditions are satisfied: 



8.2 Application 
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Table 8.1 Best (/]*) and worst (/]■ ) values of all criteria 



Criteria 





p 


IR 


IP 


OH 


UL 


OM 


tCOjjy 




Max 


Min 


Min 


Max 


Max 


Min 


Max 


ft 


56,000 


601 


1 


7,500 


25 


7.11 


9,649,680 


fr 


2,000 


5,000 


2.5 


2,000 


15 


0.60 


255,490 



1. Acceptable advantage. Q(A^) — Q{A^) > DQ, where DQ — and A^ is 
the alternative with second position in the ranking list by Q. 

2. Acceptable stability in decision-making. The alternative A^ must also be the best 
ranked by S and/or R. This compromise solution is stable within a decision making 
process, which could be the strategy of maximum group utility (when v > 0.5 is 
needed), or ‘by consensus’ (i> Rs 0.5). 

If one of the conditions is not satisfied, then a set of compromise solutions is 
proposed which consists of: 

• Alternative A' and A^ if only condition 2 is not satisfied, or 

• Alternatives A\ A^, A^ if the condition 1 is not satisfied. A^ is determined by the 
relation Q(A) — A(A^) < D Q for maximum M (the positions of these alternatives 
are in ‘closeness’). 



8.2 Application 

Let us consider as a multi-criteria problem, the same problem as in the TOPSIS 
method, the selection among 13 different RE alternatives for electric generation 
evaluated according to the same criteria (Table 7 . 1 ) . Of the criteria considered. Power, 
Operating Hours, Useful Life and Tons of emissions avoided are beneficial attributes 
and so higher values are desirable. Investment Ratio, Implementation Period and 
Operating and Maintenance Costs are non-beneficial attributes and so lower values 
are desirable. Table 8.1 shows the the best f* and the worst f.“ values of all criteria. 

The values of Sj, Rj and Qj obtained using Eqs. 8.3, 8.4, and 8.5 respectively 
are shown in Table 8.2. Ranking the alternatives, sorting by the values S, R and Q 
in decreasing order (Table 8.3) gives us, as a compromise solution, the alternative 
Ai 2 . This alternative, a Biomass plant (Co-combustion in a conventional power 
plant) of P > 50 Mw is the best ranked by the measure Q (minimum). In addition, 
conditions 1 and 2 are satisfied as this alternative is also the best ranked by S and R. 

When comparing the results obtained with the TOPSIS and VIKOR methods, 
we can see that the same alternative, the Co-combustion in a conventional power 
plant, is the highest ranked by both methods. Being the highest by the TOPSIS 
method means that this alternative is the best in terms of the ranking index. In 
addition, being the highest ranked alternative by the VIKOR method means that it 
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Table 8.2 Values of Sj, Rj and Qj 





Al 






A4 


As 


^6 


A? 


^8 


A9 


A 10 


All 


Ai2 


Ai3 


Sj 


0.713 


0.646 


0.371 


0.693 


0.621 


0.539 


0.536 


0.709 


0.698 


0.698 


0.681 


0.238 


0.545 


Rj 


0.301 


0.272 


0.183 


0.301 


0.273 


0.273 


0.266 


0.301 


0.301 


0.301 


0.301 


0.140 


0.319 


Qj 


0.955 


0.802 


0.263 


0.934 


0.778 


0.691 


0.669 


0.951 


0.939 


0.939 


0.921 


0 


0.826 



Table 8.3 Ranking 




Qj 






Rj 






Sj 




Qn 




0.000 


R\2 




0.093 


S\2 




0.131 


Qi 




0.289 


Rj 




0.151 


Sj 




0.364 


Qe 




0.330 


Re 




0.175 


Se 




0.373 


Q3 




0.582 


R3 




0.258 


S5 




0.492 


Qs 




0.611 


Rs 




0.300 


S3 




0.526 


Q 2 




0.868 


Ri 




0.383 


Sa 




0.600 


Qa 




0.882 


Ri 




0.425 


Sl 3 




0.606 


2i3 




0.922 


Ra 




0.425 


S2 




0.650 


Gii 




0.934 


Ri 




0.425 


Sn 




0.659 


G9 




0.953 


R9 




0.425 


Sio 




0.680 


Qio 




0.953 


R\o 




0.425 


S9 




0.680 


Qi 




0.963 


Rn 




0.425 


Si 




0.692 


Gs 




0.965 


R\3 




0.450 


Ss 




0.694 



is the closest to the ideal solution, and to be as close as possible to the ideal is the 
rationale of human choice. 



8.3 Concluding Remarks 



TOPSIS and VIKOR methods use different kinds of normalization to eliminate the 
units of criterion functions: the VIKOR method uses linear normalization whereas the 
TOPSIS method uses vector normalization. However, the main difference appears in 
the aggregating approaches. The VIKOR method introduces an aggregating function 
representing the distance from the ideal solution. This ranking index is an aggre- 
gation of all criteria, the relative importance of the criteria, and a balance between 
total and individual satisfaction. The TOPSIS method introduces the ranking index 
including the distances form the ideal point and from the negative-ideal point with- 
out considering their relative importance. Applications of the VIKOR method can 
be found in Sayadi et al. [4], Chang [5], Sanayei et al. [6], Bazzari et al. [7] and San 
Cristobal [8]. 
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Chapter 9 

A Multi-Criteria Data Envelopment Analysis 



To estimate the efficiency scores of companies, various methods have been developed 
during the past two decades. These methods are generally classified as parametric 
and non-parametric methods. In the parametric methods, a cost of production func- 
tion is estimated, whereas in the non-parametric methods, it is not necessary to 
estimate the cost or production function. Stochastic Frontier Analysis (SFA) and 
Data Envelopment Analysis (DBA) are the major parametric and non-parametric 
models respectively. In this chapter a Multiple Criteria Data Envelopment Analysis 
(MCDEA) model is developed and applied to the selection of a renewable project 
based on the concept of efficiency. 

9.1 DEA 

Charnes et al. [1] first introduced the DEA concept and many studies have since 
appeared that deal with various types of applications where the presence of multiple 
inputs and outputs makes comparison difficult. A non-parametric piecewise frontier 
(a best practice efficiency frontier) composed of Decision-Making Units (DMUs), 
which own the optimal efficiency over the datasets is constructed by DEA for compar- 
ative efficiency measurement. Those DMUs located at the efficiency frontier have 
their maximum outputs generated among all DMUs by taking the minimum level 
of inputs, which are efficient DMUs and own the best efficiency among all DMUs. 
The existing gap from any DMUs to the efficiency frontier shows how far the DMUs 
should be further improved to reach the optimal efficiency level. DEA produces 
detailed information on the efficiency of the unit, to be measured without any assump- 
tions regarding the functional form of the production function, not only relative to 
the efficiency frontier, but also to specihc efficient units which can be identihed as 
role models [2, 3]. Thus, DEA can be used by inefficient organizations to benchmark 
efficient and ‘best-practice organizations’. 

The classical DEA model for evaluating the efficiency of a DMU, denoted by 
DMUo is as follows [1]: 
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max /iQ — ^ryrjo 

s.t. 

zr=i>^/v/o = i (9.1) 

Zr=l >^ryrj ~ Z"Ll ViXij^ < 0 

Ur, Vi > 0 

where j is the DMU index {j = 1, n)\ r the output index (r = 1, ^); i the 
input index (i = 1, m); Xij the value of the ith input of the yth DMU\ yrj the 
value of the rth output of the y'th DMU\ Ur the weight given to the rth output; u; 
the weight given to the ith input; and /?q the relative efficiency of DMUq, the DMU 
under evaluation. In this model, DMUq is efficient if and only if ho = 1. 

Lii and Reeves [4] extend the concept of relative efficiency and the method of 
efficiency evaluation from single criterion-based conventional approach to multiple 
criteria-oriented one. The authors argue that the resulting DEA models are more 
flexible and powerful in many aspects, particularly in discriminant analysis and 
weight restriction, showing three advantages. First, in classical DEA, if a DMU is 
efficient, its optimal (weight) solution is almost surely non-unique. In this situation, 
a Linear Programming procedure returns the first optimal solution it finds usually 
giving solutions with extremely distributed weights. Second, Multi Objective Linear 
Programming (MOLP) solutions contain not only the solution that individually opti- 
mize each of the objectives, but also other non-dominated solutions that, while non- 
optimal under any given criterion, provides alternative choices for the associated 
DMU. There may be situations where some of these non-dominated solutions are 
more preferred to those solutions optimizing individual objectives. Third, the total 
number of non-dominated solutions associated with aDMf/ often reflects the stability 
of a DMUs efficiency scores relative to the changes in efficiency criteria. 

Model (9.1) can be expressed equivalently in the following deviational variable 
form: 

min do (or max ho = Xr=i “'•>'r/o) 
s.t. 

zr=i (9.2) 

X!/.= l U-rJrj — 2^! = ! ViXijQ + dj — 0 

Ur, Vi, dj > 0 

where do is the deviation variable for DMUo and dj the deviation variable for the y th 
DMU (appeared at the yth original inequality constraint). Under this model, DMUo 
is efficient if and only if r/o = 0 (or equivalently, ho = 1). If DMUo is not efficient, 
its efficiency score is /jq = \ — do. The quantity do, which is bounded by the interval 
[0, 1), can be regarded as a measure of inefficiency. The smaller the value of do, 
the less inefficient (and thus the more efficient) DMUq is. In this sense, we say that 
classical DEA method [Model (9.1) or (9.2)] minimizes DMUq’s. inefficiency, as 
measured by do, under the constraint that the weighted sum of the outputs is less 
than or equal to the weighted sum of the inputs for each DMU. 

Lii and Reeves [4] propose a MCDEA model that is different from other models 
in that each criterion is an independent objective function. That is, each criterion 
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defines a particular concept of efficiency and there is no prior preference among 
these efficiency criteria. The form of the proposed MCDEA model depends upon 
the efficiency criteria used. For a MCDEA problem that has the following three 
criteria, i.e., minimizing cLq, minimizing the sum of the deviations and minimizing 
the maximum deviation, the MCDEA model can be specified as follows: 



min r/o(or max Hq = Xr=i ^ryrjo) 
min 1:7=1 dj 
min M 
s.t. 

Z m 

1 = 1 f;Xi 7 Q = l 

Xr=l Uryrj ~ X!( = l + dj = 0 

M -dj >0 

Ur, Vi, dj > 0 



(9.3) 



The first objective of model (9.3) is identical to the objective of model (9.1) 
or (9.2). The second objective function is a straightforward representation of the 
deviation sum and the variable M in the third objective represents the maximum 
quantity among all deviation variables dj {j = !,...,«). Notice that the feasible 
region for decision variables Ur and u; in model (9.3) is the same as that in models 
(9.1) and (9.2). The effect of added constraints, M — dj > 0 is to make M the 
maximum deviation. They do not change the feasible region of decision variables. 

Model (9.3) is a Multi-objective linear programming model (MOLP) and, in a 
MOLP problem it is generally impossible to find a solution that optimizes all objec- 
tives simultaneously. Therefore, the task of a MOLP solution process is not to find 
an optimal solution but, instead, to find non-dominated solutions (in multiple criteria 
terminology, a non-dominated solution is also called an efficient solution) and to help 
select a most preferred one. One fact to point out is that a non-dominated solution set 
for a MOLP problem will always contain, but is not limited to the optimal solutions 
obtained by individually optimizing each of the objectives in the MOLP problem 
under the setting of single objective linear programming. 

The solution that optimizes the first objective function of model (9.3) is equivalent 
to the optimal solution of model (9.1) or (9.2). That is, DMUq is efficient (in the 
classical sense) if and only if the value of do corresponding to the solution that 
optimizes the first objective function of model (9.3) is zero. In the same way, it 
can be dehned a DMU’s relative efficiency corresponding to the second and third 
objective criteria in the following way: DMUq is minsum efficient and minimax 
efficient if and only if the value of do corresponding to the solution that minimizes 
the second and third objective function of model (9.3) is zero. In all the three above 
dehnitions, no matter DMUo is efficient or not, its DEA efficiency score is 1 — <fo- 

Efficiencies dehned under minsum and minimax criterion are most restrictive 
than that dehned in the classical DEA and these two criteria generally yield fewer 
efficient DMUs. That is, it is more difficult for a DMU to achieve minsum or minimax 
efficiency that to achieve classical DEA efficiency. More precisely, if DMUo is 
minsum or minimax efficient, it must also be DEA efficient, because by dehnition. 
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minsum or minimax efficiency requires do — 0. However, if DMUo is DEA efficient, 
it may or may not be minsum or minimax efficient because do does not necessarily 
imply that M or 2] dj is minimized. Including these new criteria in a DEA model 
will result in the improvement in discriminating power. On the other hand, since M 
or 2] dj are functions of all deviation variables and each deviation variable is related 
to a constraint, minimizing M or 2] dj is, in some sense, equivalent to improving 
tighter constraints on weight variables. 



9.2 Application 

In this section, the MCDEA model is applied to determine the most appropriate RE 
alternative based on the concept of efficiency [5]. The input and output data of 13 
renewable energies technologies used to perform the analysis, shown in Table 9.1 
are defined as follows: 



Inputs 


VI 


Investment ratio (Euro * 10^ /Kw) 


V2 


Implement period (years) 


V3 


Operating and maintenance costs (Euro * 10^^ /Kwh) 


Outputs 




Ml 


Power (Mw * 10^) 


U2 


Operating hours (Hours * 10^ /year) 


M 3 


Useful life (years) 


M 4 


Tons of C O 2 avoided OCO 2 * 10^ /year) 



The results of classical DEA method, solved by LP are given in Table 9.2. Using 
model (9.3), the results obtained are shown in Tables 9.3, 9.4 and 9.5, corresponding 
to the three criteria respectively (minimizing do , minimizing the sum of the deviations 
and minimizing the maximum deviation). Efficiency scores in Tables 9.2 and 9.3 are 
identical since they are obtained under the same criterion. As we can see, all DMUs 
except DMUi (Solar Thermo-electric P > lOMin) are efficient. However, we can 
easily observe differences between the solutions shown in both Tables. The input and 
output weights obtained by the MOLP procedure in Table 9.3 are distributed more 
evenly than those obtained by classical DEA in Table 9.2. 

Table 9.4 lists the efficiency scores for each DMU under the Minsum criteria. 
Under this criteria the only DMUs rated as efficient are DMU\ (Windpower P < 
5Mw), DMU 2 (Windpower 5 < P < lOMw), DMU 3 (Windpower 10 < P < 
50 Mw) and DMU\i (Biomass forest industrial wastes P < Mw). Table 9.51ists 
the efficiency scores for each DMU under the Minimax criterion. In particular, only 
DMU-i (Windpower 10 < P < 50 Mw) and DMU \2 (Biomass Co-combustion in 
conventional central P > 50 Mw) are rated as efficient, while the rest of the DMUs, 
rated as efficient by classical DEA and by Minsum criteria, are no longer efficient 
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Table 9.1 Input and output data for the 13 DMUs 



DMU 




Input weights 

Dl V2 V2 


Output weights 

U\ U2 M3 


M4 


DMUi 


Windpower P < 5 Mw 


0.937 


1 


1.47 


0.5 


2.35 


20 


1.93 


DMU2 


Windpower 5 < P < lOMiu 


0.937 


1 


1.47 


1.0 


2.35 


20 


3.22 


DMU?, 


Windpower 10 < P < 50Mw 


1.500 


1 


1.51 


2.5 


2.35 


20 


9.65 


DMUa 


Hydroelectric P < 10 Mw 


0.700 


1.5 


1.45 


0.5 


3.10 


25 


0.47 


DMUi 


Hydroelectric 10 < P < 25 Mw 


0.601 


2 


0.70 


2.0 


2.00 


25 


0.26 


DMUe 


Hydroelectric 25 < P < 50 Mw 


5.000 


2.5 


0.60 


3.5 


2.00 


25 


0.26 


DMUi 


Solar Thermo-electric P > lOMiu 


1.803 


2 


4.2 


5.0 


2.59 


25 


0.48 


DMUi 


Biomass (energetic cultivations) 
P < 5Mw 


1.803 


1 


7.11 


0.5 


7.50 


15 


2.52 


DMUg 


Biomass (forest and agricultural wastes) 
P < 5 Mill 


1.803 


1 


5.42 


0.5 


7.50 


15 


2.52 


DMUio 


Biomass (farming industrial wastes) 
P < 5 Min 


1.803 


1 


5.42 


0.5 


7.50 


15 


2.52 


DMUn 


Biomass (forest industrial wastes) 
P < 5 Mm 


1.803 


1 


2.81 


0.5 


7.50 


15 


2.52 


DMUn 


Biomass (Co-combustion in conventional 
central) P > 50 Mw 


0.856 


1 


4.56 


5.6 


7.50 


20 


4.84 


DMUi? 


Bio fuels P < 2 Mm 


1.503 


1.5 


2.51 


0.2 


7.00 


20 


5.91 



Table 9.2 Classical DEA results 


DMU 


Efficiency 


Input weights 




Output weights 






Vl 


V2 


lt3 


U\ 


»2 


M 3 


M 4 


DMUi 


1 


0.14 


0.40 


0.32 


0 


1.31 


3.46 


0 


DMU2 


1 


0.29 


0.38 


0.23 


0 


1.30 


3.44 


0 


DMU? 


1 


0 


0.44 


0.37 


0.38 


0 


0 


0 


DMUa 


1 


0.10 


0.32 


0.26 


0.01 


1.05 


2.70 


0 


DMU? 


1 


0.45 


0.31 


0.08 


0.07 


0 


3.45 


0 


DMUe 


1 


0 


0.33 


0.28 


0.28 


0 


0 


0 


DMUi 


0.78 


0 


0.20 


0.15 


0.14 


0 


0.27 


0 


DMUi 


1 


0 


1.00 


0 


0 


1.33 


0 


0 


DMUg 


1 


0 


1.00 


0 


0 


1.33 


0 


0 


DMUio 


1 


0 


1.00 


0 


0 


1.33 


0 


0 


DMUn 


1 


0 


0.12 


0.31 


0.07 


1.24 


0 


0.01 


DMUn 


1 


0 


0.20 


0.17 


0.18 


0 


0 


0 


DMUi? 


1 


0.08 


0.12 


0.28 


0.07 


1.30 


0 


0.01 



with respect to the Minimax criteria. As we can see, the only DMU rated as efficient 
under the three criteria is DMUj, (Windpower 10 < P < 50 Mw) and it can be 
considered the only non-dominated solution. 
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Table 9.3 Minimizing do results 



DMU 


Efficiency 


Input weights 




Output weights 






ft 


V2 


v? 


Ml 


M2 


M3 


M4 


DMUi 


1 


0.10 


9.02 


0 


0 


0 


5 


0 


DMU2 


1 


0 


7.89 


0.14 


0.06 


0.54 


3.64 


0 


DMU?, 


1 


0 


7.84 


0.14 


0.06 


0.54 


3.62 


0 


DMUa 


1 


0 


5.67 


0.10 


0.04 


0.39 


2.62 


0 


DMUi 


1 


0 


4.70 


0.08 


0.04 


0.32 


2.17 


0 


DMUe 


1 


0 


3.83 


0.07 


0.03 


0.26 


1.77 


0 


DMUi 


0.78 


0 


3.62 


0.07 


0.03 


0.25 


1.67 


0 


DMUi 


1 


0 


3.01 


0.10 


0.02 


0.48 


1.37 


0 


DMUg 


1 


0 


3.61 


0.12 


0.03 


0.57 


1.65 


0 


DMUw 


1 


0 


3.61 


0.12 


0.03 


0.57 


1.65 


0 


DMUn 


1 


0 


6.15 


0.12 


0.05 


0.45 


3.05 


0 


DMUn 


1 


0 


4.02 


0.13 


0.03 


0.64 


1.83 


0 


DMUn 


1 


0 


5.11 


0.09 


0.04 


0.35 


2.36 


0 



Table 9.4 Minsum DEA results 


DMU 


Efficiency 


Input weights 




Output weights 








V2 


U3 


Ml 


M2 


M3 


M4 


DMUi 


1 


0 


8.51 


0.10 


0 


0.67 


4.21 


0 


DMU2 


1 


0 


8.51 


0.10 


0 


0.67 


4.21 


0 


DMU?, 


1 


0 


8.48 


0.10 


0 


0.67 


4.19 


0 


DMUa 


0.88 


0 


5.98 


0.07 


0 


0.47 


2.96 


0 


DMUi 


0.67 


0 


4.80 


0.06 


0 


0.38 


2.37 


0 


DMUe 


0.54 


0 


3.89 


0.05 


0 


0.31 


1.92 


0 


DMUi 


0.36 


0.10 


1.10 


0.06 


0 


0.38 


1.04 


0 


DMUi 


0.53 


0.06 


1.10 


0.11 


0.02 


0.54 


0.74 


0 


DMUg 


0.65 


0.08 


1.34 


0.13 


0.03 


0.70 


0.91 


0 


DMUio 


0.65 


0 


3.61 


0.12 


0.03 


0.57 


1.65 


0 


DMUn 


1 


0.15 


3.26 


0.14 


0 


0.86 


2.36 


0 


DMUn 


0.80 


0 


1.80 


0.18 


0.02 


1.64 


1.28 


0 


DMUn 


0.86 


0 


5.56 


0.07 


0 


0.44 


2.75 


0 



9.3 Concluding Remarks 



Application of DEA in the electric power industry may be traced back to the publi- 
cations of Cote [6], Hjalmarsson and Veiderpass [7], Miliotis [8] and Golany et 
al. [9]. Since then, DEA has been widely applied both to the distribution aspects 
of operations and to the electric generation. Jha and Shrestha [10] evaluated the 
efficiency of Hydropower plants in Nepal using DEA and Vaninski [11] using DEA 
estimated the efficiency of electric power generation in the U.S. for the period of 199 1 
through 2004. Azadeh et al. [12] presented and integrated hierarchical DEA-Principal 
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Table 9.5 Minimax DEA results 



DMU 


Efficiency 


Input weights 




Output weights 






ft 


V2 


v? 


Ml 


M2 


M3 


M4 


DMUi 


0.88 


0.10 


9.02 


0 


0 


0 


5 


0 


DMU2 


0.91 


0 


7.89 


0.14 


0.06 


0.54 


3.64 


0 


DMU?, 


1 


0 


7.84 


0.14 


0.06 


0.54 


3.62 


0 


DMUa 


0.79 


0 


5.67 


0.10 


0.04 


0.39 


2.62 


0 


DMUi 


0.68 


0 


4.70 


0.08 


0.04 


0.32 


2.17 


0 


DMUe 


0.60 


0 


3.83 


0.07 


0.03 


0.26 


1.77 


0 


DMUi 


0.62 


0 


3.62 


0.07 


0.03 


0.25 


1.67 


0 


DMUi 


0.58 


0 


3.01 


0.10 


0.02 


0.48 


1.37 


0 


DMUg 


0.70 


0 


3.61 


0.12 


0.03 


0.57 


1.65 


0 


DMUw 


0.70 


0 


3.61 


0.12 


0.03 


0.57 


1.65 


0 


DMUn 


0.82 


0 


6.15 


0.12 


0.05 


0.45 


3.05 


0 


DMUn 


1 


0 


4.02 


0.13 


0.03 


0.64 


1.83 


0 


DMUn 


0.72 


0 


5.11 


0.09 


0.04 


0.35 


2.36 


0 



Component Analysis approach for location of solar plants. Sadjadi and Omrani [13] 
presented a DEA model with uncertain data for performance assessment of Iranian 
electricity distribution companies. Jayanthi et al. [14] showed a DEA-based applica- 
tion to the U.S. photovoltaic industry to evaluate the potential of innovations occur- 
ring in various stages of the industry value chain in terms of its relative efficiency 
with respect to a best practices frontier. 
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Chapter 10 

Multi- Attribute Utility Theory 



Multi-attribute utility theory is an extension of Utility theory developed to help 
decision-makers assign utility values, taking into consideration the decision-maker’s 
preferences, to outcomes by evaluating these in terms of multiple attributes and 
combining these individual assignments to obtain overall utility measures. Utility 
theory has generally been used to develop a relationship between utility and costs 
incurred as a consequence of a particular decision. There are situations where, rational 
decision-makers who subscribe the von Neuman-Morgenstern axioms are sometimes 
willing to violate the Expected Monetary Value minimization criterion (when dealing 
with benefits to maximize the Expected Monetary Value) and to sacrifice it to reduce 
risk, choosing the alternative that maximizes his or her expected utility. 



10.1 The MAU Function 

Multi-criteria utility theory generally combines the main advantages of simple scor- 
ing techniques and optimization models. Further, in situations in which satisfaction 
is uncertain, utility functions have the property that expected utility can be used as a 
guide to rational decision-making. 

Multi-attribute utility theory takes into consideration the decision-maker’s pref- 
erences in the form of utility function which is defined over a set of attributes. The 
utility value can be determined in the following three steps [1]: (i) Determination 
of single attribute utility functions. The values of utilities can vary between zero 
and one and reflect the level of importance in the achievement of that attribute; (ii) 
Verification of preferential and utility independence conditions and (iii) Derivation 
of the multi-attribute utility function. 

A utility function is a device which quantifies the preferences of a decision- 
maker by assigning a numerical index to varying levels of satisfaction of a particular 
criterion. For a single criterion, the utility of satisfaction of a consequence x’ is 
denoted by {u{x')). Utility functions are constructed such that (m(jc')) is less preferred 
to (u(x")), i. e. 
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Fig. 10.1 Utility functions 



Risk aversion 




u{x') < u(x") 

if and only if x' is less preferred to x”, i. e. x' < x" 

In other words, a utility function is a transformation of some level of project 
(alternative) performance, x' , measured in its natural units into an equivalent level of 
decision-maker satisfaction. Utility functions contain information about the decision- 
maker’s attitude toward risk being reflected in the shape of the utility curve which 
combines the decision-maker’s preference attitudes, i.e. increasing or decreasing 
utility with increasing or decreasing x’ . Depending on the decision maker’s attitudes 
toward risk, utility functions can be concave, convex or linear as shown in Fig. 10.1. 
Concave utility functions describe risk-averse situations, convex utility functions 
describe risk- seeking situations and linear utility functions describe a risk-neutral 
situation. 

All decisions involve choosing one, from several alternatives. Typically, each 
alternative is assessed for desirability on a number of scored criteria. What con- 
nects the criteria scores with desirability is the utility function. The most common 
formulation of a utility function is the additive model; 



{/, = Wj, f/,y for alii (10.1) 

where t/, is the overall utility value of alternative /; Uij is the utility value of the yth 
criterion for the ith alternative and Uij equals u(Xi), for 1 > i > n and 1 > j > m. 
Xi designates a specihc value of x,y ; n is the total number of criteria, m is the total 
number of alternatives and W j is the relative weight of the yth criterion. 

The advantage of the additive form is its simplicity. In order to determine the 
overall utility function for any alternative, a decision-maker needs to determine only 
the n utility functions for that alternative. For the construction of the utility functions 
the decision-maker’s preferences are analyzed by the method suggested by Bell et al. 
[2], and Keeney and Raiffa [3]. The hrst step involves the identification of the best 
and the worst outcomes (criteria scores) for each one of the criteria. The decision- 
maker is free to set these utility values at any level provided that the best outcome has 
the higher value. The usual method is to assign the worst outcome a utility value of 
zero and the best outcome a utility value of unity. This establishes the range of utility 
values to from 0 to 1 between the worst and the best possible outcomes. To determine 
the utility of intermediate values, the decision-maker is offered the following options: 
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p = ? 




Rl 



18 points U = 1 



10 points U = 0 



13 points U = ? 
16 points 
15 points 
11 points 



Fig. 10.2 Routes 



1 . Certain option. In this case the decision-maker is offered a certain outcome with 
a probability p = 1 ; 

2. Risk option. In this case the decision-maker is offered a probabilistic outcome 
in the form of a gamble, in which the decision-maker either receives the best 
outcome with a probability p or the worst outcome with a probability of 1 — p. 

The following is an example of how the utility values for an hypothetical criterion 
are obtained and from which the utility curve is established. Let us consider the 
following scores of five alternatives for an hypothetical criterion: 





A 


B 


C 


D 


E 


Scores 


13 


14 


10 


18 


16 



The first step is to identify the best and worst outcomes for the criterion considered 
and assign arbitrary utility values of 1 for the best outcome (alternative D with 18 
points) and 0 for the worst outcome (alternative C with 10 points). The utility of 
the intermediate values is then determined by offering the decision-maker a choice 
between the following lotteries (Fig. 10.2): 

1. Lottery 1: go to route Ri for a certain consequence of 13 points (alternative A). 

2. Lottery 2: go to route R 2 for either a best consequence of 18 points (alternative 
D) with a probability of p or a worst consequence of 10 points (alternative C) 
with a probability of 1 — p. 

What utility value should the decision-maker assign to a certain outcome of 13 
points? For the decision-maker to make good decision and choose from the two 
routes, the utility value of 13-point score must be assessed and compared with the 
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Table 10.1 Significant criteria in cooking energy decisions 



Number 


Criteria 


Desired value 


Technical 






Cl 


Fuel consumption 


Low 


C2 


Cooking time 


Low 


C3 


Durability 


High 


C4 


Quality, reliability 


High 


Cs 


Sophistication level 


High 


Ce 


Size/Weight/Space needs 


Low 


Cl 


Ruggedness 


High 


Cg 


Continuity of use 


High 


C9 


Need for tracking 


Low 


Cio 


Nutrition value of food 


High 


Economic 






Cii 


Initial cost 


Low 


Cl 2 


Fuel cost per month 


Low 


Cn 


Maintenance cost 


Low 


Cl 4 


Available subsidy 


Low 


Cl 5 


Rate of interest on loan, if any 


Low 


Social 






C16 


Pollution hazards 


Low 


Cn 


Fluman drudgery 


Low 


Cis 


Overall safety 


High 


Behavioral 






C19 


Aesthetics 


High 


C20 


Motivation to buy 


High 


C2I 


Taste of food 


High 


C22 


Cleanliness of utensils 


High 


C23 


Ease of operation 


High 


C24 


Type of dishes cooked 


High 


C25 


Need for additional cooking system 


Low 


Commercial 






C26 


Improvements in models 


High 


C27 


Spares and after sales services 


High 


C28 


Distribution network 


High 


C29 


Market research 


High 


C30 


Need for user training 


Low 



expected utility of the risk option. To do this, the decision-maker determines a rela- 
tive preference for a 13-point consequence by finding the probability p for the best 
outcome, to which the decision-maker is indifferent, between the certain route R\ for 
a 13-point outcome and the gamble route R 2 for the two possible outcomes of 18 and 
10 points. Let us assume that there is a probability of 0.3 for getting the best outcome 
and a probability of 0.7 for getting the worst outcome from the route R 2 ■ Which route 
would the decision-maker prefer in this case? Since/? = 0.3 the chance of getting the 
best outcome from route R 2 is very small, so in this case the decision-maker will 
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Table 10.2 Alternatives 



Alternative 


Description 


At 


Chulha 


A2 


Improved chulha 


^3 


Kerosene stove 


A4 


Biogas stove 


45 


LPG stove 


^6 


Microwave oven 


A7 


Electric oven 


^8 


Solar box cooker 


A9 


Parabolic solar cooker 



not gamble. He prefers to choose route Ri with a 13-point certain outcome. Now, 
let us assume that there is a probability of 0.9 for getting the best outcome and a 
probability 0.1 for getting the worst outcome from route /? 2 - Which route does the 
decision-maker now prefer?. Since p = 0.9, in this case there is a high chance of get- 
ting the best outcome of 18 points, so he will gamble and choose route Rj. Now, 
let us take a probability of 0.45 for getting the best outcome and a probability of 
0.55 for getting the worst outcome from route R 2 - Which route does he or she now 
prefer? Putting p = 0A5 makes the thing difficult to choose for the decision-maker, 
but he will go for the certain outcome route R^. Doing some more of these trials and 
errors, the decision-maker considers that a probability of 0.5 will make him indiffer- 
ent between the two routes R \ and Rj ■ According to utility theory, by choosing the 
probability that makes him indifferent between the two routes, the decision-maker 
has assigned a utility value for the certain outcome of 13 points. It is known from 
the principles of probabilities that the expected value of any random variable in the 
space will equal the sum of probability of each variable times its score. In this case, 
the expected utility for the route R 2 which includes two variables or two outcomes 
(the best outcome with m = 1 and the worst outcome with m = 0) will be: 

p (utility of best outcome) -|- (1 — /?) (utility of worst outcome) 

= 0.5 * t/(18) -t- (1 - 0.5)17(10) = 0.5 * 1 -b 0.5 * 0 

Since the decision-maker is indifferent between 13 points for certain and this 
gamble, the alternatives must have the same utility value that is, t/(13) = 0.5. This 
procedure can be used for any scores between 10 and 18. The more utility values 
obtained, the better the utility curves appear. 



10.2 Application 



Let us consider the case by Pohekar and Ramachandran [ 1 ] of a parabolic solar cooker 
(PSC) which is evaluated with respect to eight prevalent domestic cooking devices in 
India. Thirty different criteria categorized under technical, economic, environmental. 
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Table 10.3 Evaluation matrix for all criteria 



Criteria 


Ai 


^2 


M 


^4 


As 




Ai 


As 


A9 


Scaling constant 


Cl 


2 


1 


0.5 


0 


0.25 


2 


2 


0 


0 


0.0355 


Cl 


60 


60 


30 


15 


15 


5 


30 


180 


20 


0.0384 


Cl 


1 


4 


15 


5 


20 


5 


5 


10 


20 


0.0397 


C 4 


2 


2 


8 


4 


10 


10 


6 


6 


6 


0.0376 


Cs 


2 


2 


6 


6 


10 


10 


10 


4 


4 


0.0326 


Ce 


2 


1 


2 


50 


10 


5 


3 


5 


15 


0.0332 


Cl 


2 


2 


6 


6 


10 


8 


8 


8 


8 


0.0310 


Cs 


10 


10 


10 


4 


10 


10 


10 


4 


2 


0.0335 


C 9 


0 


0 


0 


0 


0 


0 


0 


4 


10 


0.0326 


Cio 


6 


6 


6 


6 


6 


2 


2 


10 


8 


0.0371 


Cii 


10 


50 


200 


5000 


4000 


8000 


5000 


2000 


7000 


0.0371 


C \2 


20 


10 


100 


0 


250 


200 


400 


0 


0 


0.0293 


Cl 3 


0 


0 


50 


200 


50 


200 


200 


50 


20 


0.0274 


Cl 4 


0 


50 


0 


2000 


0 


0 


0 


500 


2000 


0.0252 


C\s 


0 


0 


0 


0 


13 


13 


13 


3 


3 


0.0230 


Cl6 


10 


10 


8 


2 


4 


10 


10 


0 


0 


0.0217 


C\i 


10 


10 


6 


6 


2 


2 


2 


2 


2 


0.0296 


Cis 


2 


2 


4 


10 


6 


2 


2 


10 


8 


0.0433 


Cl 9 


2 


2 


4 


4 


10 


10 


8 


8 


8 


0.0348 


Cio 


2 


2 


6 


6 


10 


8 


4 


4 


4 


0.0314 


C21 


6 


6 


2 


6 


6 


2 


10 


10 


10 


0.0381 


C22 


2 


2 


2 


8 


10 


10 


10 


10 


10 


0.0368 


C23 


10 


10 


10 


4 


8 


4 


4 


6 


4 


0.0387 


C24 


10 


10 


10 


10 


10 


2 


2 


2 


4 


0.0332 


C25 


2 


2 


4 


4 


6 


8 


10 


10 


10 


0.0306 


C26 


2 


2 


6 


4 


10 


10 


6 


6 


6 


0.0342 


Cn 


10 


10 


10 


4 


10 


6 


4 


4 


2 


0.0352 


C28 


10 


6 


6 


6 


2 


6 


2 


2 


2 


0.0355 


C29 


2 


2 


8 


2 


10 


10 


4 


4 


4 


0.0315 


C30 


2 


2 


4 


8 


4 


10 


6 


6 


8 


0.0323 



social behavioral and commercial aspects are considered for the evaluation. To illus- 
trate the use of different utility functions, the model presented by the authors, which 
only uses linear utility functions, is extended to concave and convex utility func- 
tions showing decision-maker’s risk-aversion or risk-seeking attitude toward risk. 
Table 10. 1 presents a summary of 30 criteria (twenty one criteria of qualitative nature 
and nine of qualitative nature) and their significance in evaluation of a cooking device. 
Table 10.2 shows the alternatives to evaluate and Table 10.3 the evaluation matrix 
for all criteria. The values of 2, 4, 6, 8 and 10 for qualitative criteria are allotted 
to indicate values of performance of a device as very low, low, moderate, high and 
very high, respectively. Scaling constants indicate weighted average importance of 
criteria. 
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Table 10.4 Normalized evaluation matrix linear utility function, ri.sk neutral 



Criteria 


Ai 


^2 


A3 


A4 


As 


Ae 


A7 


^8 


A9 


Cl 


0.000 


0.500 


0.750 


1.000 


0.850 


0.000 


0.000 


1.000 


1.000 


Cl 


0.686 


0.686 


0.857 


0.0943 


0.943 


1.000 


0.857 


0.000 


0.914 


Cj 


0.000 


0.158 


0.737 


0.210 


1.000 


0.210 


0.210 


0.474 


1.000 


C4 


0.000 


0.000 


0.750 


0.250 


1.000 


1.000 


0.500 


0.500 


0.500 


C5 


0.000 


0.000 


0.500 


0.500 


1.000 


1.000 


1.000 


0.250 


0.250 


Ce 


0.980 


1.000 


0.980 


0.000 


0.816 


0.918 


0.960 


0.918 


0.714 


Cl 


0.000 


0.000 


0.500 


0.500 


1.000 


0.750 


0.750 


0.750 


0.750 


Cs 


1.000 


1.000 


1.000 


0.250 


1.000 


1.000 


1.000 


0.250 


0.000 


C9 


1.000 


1.000 


1.000 


1.000 


1.000 


1.000 


1.000 


0.600 


0.000 


Cio 


0.500 


0.500 


0.500 


0.500 


0.500 


0.000 


0.000 


1.000 


0.750 


Cii 


1.000 


0.995 


0.976 


0.375 


0.501 


0.000 


0.375 


0.751 


0.125 


Cl2 


0.950 


0.975 


0.750 


1.000 


0.375 


0.500 


0.000 


1.000 


1.000 


Cl3 


1.000 


1.000 


0.750 


0.000 


0.750 


0.000 


0.000 


0.750 


0.900 


Cl4 


1.000 


0.980 


1.000 


1.000 


0.000 


0.000 


0.000 


0.750 


1.000 


Cl5 


1.000 


1.000 


1.000 


0.769 


0.000 


0.000 


0.000 


0.769 


0.769 


Cl6 


0.000 


0.000 


0.200 


0.800 


0.600 


0.000 


0.000 


1.000 


1.000 


Cl7 


0.000 


0.000 


0.500 


0.500 


1.000 


1.000 


1.000 


1.000 


1.000 


Cis 


0.000 


0.000 


0.250 


1.000 


0.500 


0.000 


0.000 


1.000 


0.750 


Cl9 


0.000 


0.000 


0.250 


0.250 


1.000 


1.000 


0.750 


0.750 


0.750 


Cio 


0.000 


0.000 


0.500 


0.500 


1.000 


0.750 


0.250 


0.250 


0.250 


C2I 


0.500 


0.500 


0.000 


0.500 


0.500 


0.000 


1.000 


1.000 


1.000 


C22 


0.000 


0.000 


0.000 


0.750 


1.000 


1.000 


1.000 


1.000 


1.000 


C23 


1.000 


1.000 


1.000 


0.000 


0.750 


0.000 


0.000 


0.330 


0.000 


C24 


1.000 


1.000 


1.000 


1.000 


1.000 


0.000 


0.000 


0.000 


0.250 


C25 


1.000 


1.000 


0.750 


0.750 


0.500 


0.250 


0.000 


0.000 


0.000 


C26 


0.000 


0.000 


0.500 


0.250 


1.000 


1.000 


0.500 


0.500 


0.500 


Cn 


1.000 


1.000 


1.000 


0.250 


1.000 


0.500 


0.250 


0.250 


0.000 


C28 


1.000 


0.500 


0.500 


0.500 


0.000 


0.500 


0.000 


0.000 


0.000 


C29 


0.000 


0.000 


0.750 


0.000 


1.000 


1.000 


0.250 


0.250 


0.250 


C30 


1.000 


1.000 


0.750 


0.250 


0.750 


0.000 


0.500 


0.500 


0.750 



Table 10.4 shows the normalized evaluation matrix and utility for all the criteria 
and the alternatives considered using linear utility functions. The best outcome of 
each one of the criteria has been assigned a utility value of 1 and the worst outcome 
a utility value of zero. Since there are two kinds of criteria; the maximization criteria 
(the maximum value is desirable) and the minimization criteria (the minimum value 
is desirable) the intermediate utility values are obtained by normalizing the evaluation 
matrix as follows: 

1 . For maximization criteria 



"^min 
Amax - 



( 10 . 2 ) 
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Table 10.5 Normalized evaluation matrix concave utility function, risk aversion 



Criteria 


Ai 


^2 


A3 


A4 


As 


Ae 


A? 


M 


A9 


Cl 


0.249 


0.944 


1.000 


0.929 


0.999 


0.249 


0.249 


0.929 


0.929 


Cl 


1.000 


1.000 


1.000 


0.967 


0.967 


0.929 


1.000 


0.249 


0.982 


Cl 


0.249 


0.545 


1.000 


0.628 


0.929 


0.628 


0.628 


0.925 


0.929 


C4 


0.249 


0.249 


1.000 


0.685 


0.929 


0.929 


0.944 


0.944 


0.944 


Cs 


0.249 


0.249 


0.944 


0.944 


0.929 


0.929 


0.929 


0.685 


0.685 


Ce 


0.943 


0.929 


0.944 


0.249 


1.000 


0.980 


0.957 


0.980 


1.000 


Cl 


0.249 


0.249 


0.944 


0.944 


0.929 


1.000 


1.000 


1.000 


1.000 


Cs 


0.929 


0.929 


0.929 


0.685 


0.929 


0.929 


0.929 


0.685 


0.249 


C9 


0.929 


0.929 


0.929 


0.929 


0.929 


0.929 


0.929 


0.998 


0.249 


Cio 


0.944 


0.944 


0.944 


0.944 


0.944 


0.249 


0.249 


0.929 


1.000 


Cii 


0.929 


0.933 


0.946 


0.837 


0.944 


0.249 


0.837 


1.000 


0.490 


C\2 


0.962 


0.947 


1.000 


0.929 


0.837 


0.944 


0.249 


0.929 


0.929 


Cl3 


0.929 


0.929 


1.000 


0.249 


1.000 


0.249 


0.249 


1.000 


0.989 


Cl4 


0.929 


0.943 


0.929 


0.929 


0.249 


0.249 


0.249 


1.000 


0.929 


C\s 


0.929 


0.929 


0.929 


1.000 


0.249 


0.249 


0.249 


1.000 


1.000 


Cl6 


0.249 


0.249 


0.612 


1.000 


0.998 


0.249 


0.249 


0.929 


0.929 


C\i 


0.249 


0.249 


0.944 


0.944 


0.929 


0.929 


0.929 


0.929 


0.929 


Cis 


0.249 


0.249 


0.685 


0.929 


0.944 


0.249 


0.249 


0.929 


1.000 


Cl9 


0.249 


0.249 


0.685 


0.685 


0.929 


0.929 


1.000 


1.000 


1.000 


Cio 


0.249 


0.249 


0.944 


0.944 


0.929 


1.000 


0.685 


0.685 


0.685 


C 21 


0.944 


0.944 


0.249 


0.944 


0.944 


0.249 


0.929 


0.929 


0.929 


C22 


0.249 


0.249 


0.249 


1.000 


0.929 


0.929 


0.929 


0.929 


0.929 


C23 


0.929 


0.929 


0.929 


0.249 


1.000 


0.249 


0.249 


0.787 


0.249 


C24 


0.929 


0.929 


0.929 


0.929 


0.929 


0.249 


0.249 


0.249 


0.685 


C 25 


0.929 


0.929 


1.000 


1.000 


0.944 


0.685 


0.249 


0.249 


0.249 


C26 


0.249 


0.249 


0.944 


0.685 


0.929 


0.929 


0.944 


0.944 


0.944 


Cn 


0.929 


0.929 


0.929 


0.685 


0.929 


0.944 


0.685 


0.685 


0.249 


C28 


0.929 


0.944 


0.944 


0.944 


0.249 


0.944 


0.249 


0.249 


0.249 


C29 


0.249 


0.249 


1.000 


0.249 


0.929 


0.929 


0.685 


0.685 


0.685 


C30 


0.929 


0.929 


1.000 


0.685 


1.000 


0.249 


0.944 


0.944 


1.000 



1 . For minimization criteria 



^max — Aj 
Amax - 



(10.3) 



where A / represents the score assigned to cooking device in evaluation matrix. 
Tlmax and Ajj^jjj are the maximum and minimum scores assigned for the selected 
criteria for the identified devices. 

Tables 10.5 and 10.6 show utility values using concave (exhibiting risk-averse 
behavior) and convex utility functions (exhibiting risk-seeking behavior), respec- 
tively. The final utilities of the devices are computed by multiplying the normalized 
utility values by respective scaling constants and the final utility value is found by 
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Table 10.6 Normalized evaluation matrix convex utility function, risk seeking 



Criteria 


Ai 


^2 


A 3 


A 4 


As 


Ae 


A? 


^8 


A 9 


Cl 


0.050 


0.138 


0.407 


0.827 


0.598 


0.050 


0.050 


0.827 


0.827 


Cl 


0.324 


0.324 


0.569 


0.718 


0.718 


0.827 


0.569 


0.050 


0.666 


Cl 


0.050 


0.013 


0.389 


0.014 


0.827 


0.014 


0.014 


0.118 


0.827 


C4 


0.050 


0.050 


0.407 


0.019 


0.827 


0.827 


0.138 


0.138 


0.138 


Cs 


0.050 


0.050 


0.138 


0.138 


0.827 


0.827 


0.827 


0.019 


0.019 


Ce 


0.19,1 


0.827 


0.787 


0.050 


0.504 


0.674 


0.748 


0.674 


0.360 


Cl 


0.050 


0.050 


0.138 


0.138 


0.827 


0.407 


0.407 


0.407 


0.407 


Cs 


0.827 


0.827 


0.827 


0.019 


0.827 


0.827 


0.827 


0.019 


0.050 


C9 


0.827 


0.827 


0.827 


0.827 


0.827 


0.827 


0.827 


0.228 


0.050 


Cio 


0.138 


0.138 


0.138 


0.138 


0.138 


0.050 


0.050 


0.827 


0.407 


Cii 


0.827 


0.817 


0.781 


0.060 


0.138 


0.050 


0.060 


0.409 


0.015 


C\2 


0.731 


0.779 


0.407 


0.827 


0.059 


0.138 


0.050 


0.827 


0.827 


Cl3 


0.827 


0.827 


0.407 


0.050 


0.407 


0.050 


0.050 


0.407 


0.641 


Cl4 


0.827 


0.788 


0.827 


0.827 


0.050 


0.050 


0.050 


0.407 


0.827 


C\s 


0.827 


0.827 


0.827 


0.434 


0.050 


0.050 


0.050 


0.434 


0.434 


Cl6 


0.050 


0.050 


0.013 


0.479 


0.228 


0.050 


0.050 


0.827 


0.827 


C\i 


0.050 


0.050 


0.138 


0.138 


0.827 


0.827 


0.827 


0.827 


0.827 


Cis 


0.050 


0.050 


0.019 


0.827 


0.138 


0.050 


0.050 


0.827 


0.407 


Cl9 


0.050 


0.050 


0.019 


0.019 


0.827 


0.827 


0.407 


0.407 


0.407 


Cio 


0.050 


0.050 


0.138 


0.138 


0.827 


0.407 


0.019 


0.019 


0.019 


C 21 


0.138 


0.138 


0.050 


0.138 


0.138 


0.050 


0.827 


0.827 


0.827 


C 22 


0.050 


0.050 


0.050 


0.407 


0.827 


0.827 


0.827 


0.827 


0.827 


C 23 


0.827 


0.827 


0.827 


0.050 


0.407 


0.050 


0.050 


0.040 


0.050 


C 24 


0.827 


0.827 


0.827 


0.827 


0.827 


0.050 


0.050 


0.050 


0.019 


C 25 


0.827 


0.827 


0.407 


0.407 


0.138 


0.019 


0.050 


0.050 


0.050 


C 26 


0.050 


0.050 


0.138 


0.019 


0.827 


0.827 


0.138 


0.138 


0.138 


Cn 


0.827 


0.827 


0.827 


0.019 


0.827 


0.138 


0.019 


0.019 


0.050 


C 28 


0.827 


0.138 


0.138 


0.138 


0.050 


0.138 


0.050 


0.050 


0.050 


C 29 


0.050 


0.050 


0.407 


0.050 


0.827 


0.827 


0.019 


0.019 


0.019 


C 30 


0.827 


0.827 


0.407 


0.019 


0.407 


0.050 


0.138 


0.138 


0.407 



Table 10.7 Results 


Utility 


Ai 


^2 


^3 


A 4 


As 


^6 


A? 


As 


A 9 


Risk neutral 


0.508 


0.515 


0.660 


0.513 


0.759 


0.485 


0.416 


0.580 


0.563 


Risk aversion 


0.635 


0.672 


0.878 


0.787 


0.887 


0.650 


0.641 


0.808 


0.762 


Risk seeking 


0.410 


0.388 


0.404 


0.286 


0.534 


0.365 


0.279 


0.358 


0.373 



adding the utilities of the devices on different criteria as follows: 

30 

Utility = ^ W jUij 
1 



( 10 . 4 ) 
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The results, ranking the alternatives considered, are shown in Table 10.7. The 
utility of the LPG stove alternatives (A 5 ) is found to be the highest, followed by 
the Kerosene stove (A 3 ), for the risk-neutral and risk-aversion situations and by the 
Chulha (Ai) for the risk-seeking situation. 



10.3 Concluding Remarks 

Multi-attribute utility theory provides one such approach and is specially useful as it 
allows the treatment of both quantitative and qualitative criteria. Applications of this 
method can be found in the selection of solar photovoltaic alternatives [4], energy 
policy making [5], environmental impact assessment [ 6 ] and electric power system 
expansion planning [7]. A more detailed discussion of Multi-attribute utility theory 
can be found in Bunn [ 8 ] and Keeney and Raiffa [9]. 
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Chapter 11 

Fuzzy PROMETHEE 



In many decision-making problems the decision-maker’s judgments are not crisp, 
and it is relatively difficult for the decision-maker to provide precise numerical val- 
ues for the criteria or attributes. These kinds of criteria make the evaluation process 
hard and vague. To deal with vagueness of human thought Zadeh [1] first intro- 
duced the fuzzy set theory, which was oriented to the rationality of uncertainty due 
to imprecision or vagueness. A major contribution of fuzzy set theory is its capa- 
bility of representing vague data. In a classical set, an element belongs to, or does 
not belong to, a set whereas an element of a fuzzy set naturally belongs to the set 
with a membership value from the interval [0,1]. In this chapter, the procedure of 
the PROMETHEE method described in Chap. 5 will he applied making the assump- 
tion that the performance of alternative solutions are fuzzy while the preferences of 
the decision-maker, such as the parameters of generalized criteria and the weighting 
factors, are not. That is, the performance of alternative solutions can be determined 
only approximately and therefore is introduced into the calculations as a fuzzy num- 
ber. 



11.1 Fuzzy PROMETHEE 

Fuzzy numbers can be presented, according to Dubois and Prade [2], in the form: 
X = (m, a,b)iR. When the variable x has the value m it is certain to belong to the 
specific class and its membership function is fix) — 1. If it has values smaller than 
(m — a) and bigger than (m + b) it does not belong in the specific class and /(x) = 0. 
In the interval between: m — a < x < m + b there is a ‘grade of membership’ of x in 
the specific class associated with a number between 0 and 1 . L and R are the functions 
indicating the change of /(x) withx to the left and right of m, respectively. Assuming 
linear functions for L and R, the fuzzy number x = (1, 0.1, 0.1) is presented in 
Fig. 11.1 as a triangle. The number 1.05 has a membership grade of 0.5 in the set. 
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Fig. 11.1 Fuzzy number X = (1,0.1, 0.1) 



Table 11.1 Operations with fuzzy numbers (m, a, b)iR 



Addition 
Opposite 
Substraction 
Multiplication by scalar 
Multiplication by fuzzy: 

- For m > 0, n > 0 

- For m < 0, n > 0 

- For m < 0, n < 0 
Inverse (for m > 0) 



(m, a, b)iR © (n, c, d)iR = (m + n, a + c, b + d)iR 
-(m, a, b)LR = (~m, b, o)rr 

{m, a, b)iR - {n, c, d)Ri = (m - n,a +d, b + c)rr 
{m, a, b)iR X {n, 0, 0)lr = (mn, an, bn)iR 



(m, a, b)iR ® (n, c, d)iR 
(m, a, b)RL ® (n, c, d)iR ' 
(m, a, b)iR ® (n, c, d)iR : 



(mn, cm + an, dm + bn)iR 
(mn, an — dm, bn — cm)Ri 
(mn, —bn — dm, —an — cm)Ri 



(m,a,b)j\,^ (m ^ , bm ^,am ^)rl 



The basic operations with fuzzy numbers used here are summarized in Table 11.1. 



11.2 Application 

Let us consider multi-criteria problem consisting in the selection of four alterna- 
tive exploitation alternatives for the low enthalpy geothermal field presented in 
Chap. 5 [3]. Considering that the relevant data cannot be determined accurately, 
for the evaluation by the Fuzzy PROMETHEE method, the following statements are 
expressed as fuzzy inputs: 

1 . Economic evaluation. The net present value can deviate from the estimated value 
in the range of ±10 percent for alternative 1; —15 percent to ±10 percent for 
alternative 2; —18 percent to ±10 percent for alternative 3 and —20 percent to 
±10 percent for alternative 4. 

2. Jobs. The number of job positions can deviate by ±2 for alternative 1 and by ±4 
for alternatives 2, 3 and 4. 

3. Energy consumption. It can deviate by ±20 percent from the estimated value. 

4. Risk index. In the range of — 1 percent to ±4 . 
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Table 11.2 Performance of alternatives expressed as fuzzy numbers [3] 



Scenarios 


NPV 


Jobs 


Energy use 


Risk 


Ai 


(373,37,37) 


(44,2,2) 


(22.4,4.48,4.48) 


(3,1.2) 




(3706,556,371) 


(99,4,4) 


(51.0,10.2,10.2) 


( 6 , 1 . 2 ) 


Ai 


(3809,685,381)) 


(100,4,4) 


(49.3,9,86,9,86) 


(8,1,3) 


M 


(3860,772,386) 


(104,4,4) 


(47.5,9.5,9.5) 


(9,1,4) 



Table 11.3 Values of 
n(Ai,Ak) 



Table 11.4 Data for 
preorders P+ and P 





Ai 




^3 


A 4 


Ai 


- 


0.375 


0.500 


0.500 


A 2 


0.500 


- 


0.125 


0.125 


^3 


0.500 


0.825 


- 


0 


A 4 


0.500 


0.457 


0.315 


- 







Ai 


^2 


^3 


A 4 


<P+(Ai) 

<p-(Ak) 


1.375 

1.500 


0.750 

0.915 


0.707 

0.940 


1.272 

0.750 



According to the above estimations, the performance of the four alternatives 
expressed as fuzzy numbers is shown in Table 1 1 .2. 

To compare these fuzzy numbers the index proposed by Yager [4] will be used. 
The index corresponds to a ‘weight average’ of the fuzzy number and is determined 
by the center of weight of the surface representing its membership function; 



F(m, a, b) 



(3m — a + b) 
3 



( 11 . 1 ) 



Following step by step the same procedure of the PROMETHEE method described 
in Chap. 5 for the case where the performance of alternative solutions are considered 
with precise numerical values, we obtain for every couple of actions the values 
jt(Ai, Ak) shown in Table 11.3 and in Table 11.4 the outgoing and the incoming 
flow. 

The intersection of the preorders is: 



A4P<‘’A3, A2P^‘U3 



Supposing that the decision-maker requests a total preorder, the net flows are 
calculated. Alternative A 4 is, again, the best ranked alternative. However, mak- 
ing the assumption that the performance of alternative solutions is fuzzy has 
changed the total preorder with regard to the PROMETHEE method as shown in 
Table 11.5. In this example criteria Jobs and Energy consumption have the same 
contribution to the outcome whether scalar numbers or fuzzy numbers are used 
because the Yager index of the fuzzy number coincides with m. The difference in 
performance arises from the two other criteria. Economic evaluation and Risk. The 
combined differences in these two criteria resulted in a change in the total ranking 
of the alternatives. 
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Table 11.5 Ranking 



PROMETHEE A4 > A3 > Ai > A2 

F-PROMETHEE A4 > Ai > A2 > A3 



11.3 Concluding remarks 

The Fuzzy PROMETHEE method is preferred when substantial uncertainties and 
subjectivities exist in certain environments and the vagueness of the linguistic terms 
in the evaluation process requires the employment of fuzzy numbers. In the literature, 
the method has been applied in environmental management for ranking contaminated 
sites [5], in logistic service selection [6] and material handling equipment selection 
[7, 8]. In the renewable energy industry, the applications are limited to the work by 
Goumas and Lygerou [3]. 
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Chapter 12 

FAHP 



The AHP method introduced by Saaty [1] shows the process of determining the 
priority of a set of alternatives and the relative importance of attributes in a multi- 
criteria decision-making problem. The primary advantage of the AHP approach is 
the relative ease with which it handles multiple criteria and performs qualitative and 
quantitative data. However, AHP is frequently criticized for its inability to adequately 
accommodate the inherent uncertainty and imprecision associated with mapping 
decision-maker perceptions to extract numbers [2] . 

Analysis of hierarchical structures in fuzzy environment was initially proposed 
by Buckley [3], who examined expressions of decision makers regarding with the 
pairwise comparisons while utilizing fuzzy ratios instead of crisp values. There are 
many alternatives as solution approaches [4, 5, 6] to perform on the Fuzzy AHP- 
based study on multiple criteria decision making problems. Buyukozkan et al. [7] 
introduced the main characteristics of different approaches on solving FAHP-based 
models by expressing their advantages and disadvantages with several points of view. 
In this chapter we use two methods in the selection of RE projects. The first one is a 
modified analytical hierarchy process proposed by Zeng et al. [8] while the second 
one is based on Chang’s extent analysis [9, 10]. 



12.1 FAHP 

The FAHP methodology consists of the following steps [8]: 

Step 1. Measure factors in the hierarchy 

The experts are required to provide their judgments on the basis of their knowledge 
and expertise for each factor at the bottom level in the hierarchy. Each expert can 
provide his or her judgments as a precise numerical value, a range of numerical 
values, a linguistic term or a fuzzy number. A score system is shown in Fig 12.1. 
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Fig. 12.1 Score system 



U(x) 




Fig. 12.2 Score system 




Step 2. Compare factors using pair-wise comparisons 

The experts are required to compare every factor pair-wise in their corresponding 
section structured in the hierarchy and calibrate them on either a crisp or a fuzzy 
scale. 

Step 3. Convert preferences into Standardized Trapezoidal 
Fuzzy Numbers (STFNs) 

Because the values provided by experts are crisps, e.g., a numerical value, a range 
of numerical value, a linguistic term or a fuzzy number, the STFN is employed to 
convert these experts’ judgments into a universal format for the composition of group 
preferences. A STFN can be defined as A = (a, b, c, d), where Q<a<b<c<d 
as shown in Fig. 12.2, and its membership function is defined as follows: 

fora < X < b 

(b-a) — — 

1 forb < X < c 

( 12 . 1 ) 

forc<x<d 
0 for otherwise 
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Step 4. Aggregate individual STFNs into group STFNs 

The aim of this step is to apply an appropriate operator to aggregate individual 
preferences made by individual experts into a group preference of each factor. The 
aggregation of STFN scores is performed by applying the fuzzy weighted trapezoidal 
averaging operator, which is defined by: 

S — Sil ® Cl © 5/2 ® C 2 © • • • © Sim ® c,„ (12.2) 

where 5/ is the fuzzy aggregated score of the factor Fj, Sn, 5/2, • ■ . , 5/m are the 
STFN scores of the factor f, measured by m experts E\, E 2 , . ■ ■ , Em respectively, 
® and © denote the fuzzy multiplication operator and the fuzzy addition operator 
respectively, and ci, C 2 , . . . , Cm are contribution factors (CFs) allocated to experts, 
El, E 2 , ■ ■ ■ , Em and ci + C 2 + • • ■ + c,„ = 1. Similarly, the aggregation of STFN 
scales is defined as 



a — an ® ci © a /2 <S> C 2 © • • • © a/m ® c,„ (12.3) 

wherea/j is the aggregated fuzzy scale of F/ comparing to F,- ; i, j — \,2, . . ., n; a/ji, 
aij 2 , . . ■ , a/ym are the corresponding STFN scales of F, comparing to Fj measured 
by experts E\, E 2 , ■ ■ . , Em, respectively. 

Step 5. Defuzzify the SFTN scales 

In order to convert the aggregated STFN scales into matching crisp values that 
can adequately represent the group preferences, a proper defuzzification method is 
needed. Assuming an aggregated STFN scale atj — {a\j, a™, a"j, a^S), the matching 
crisp value aij can be obtained 



an 



a|.+2(a;;;+a«) + a“ 



(12.4) 



where a/y = 1, aji — Consequently, all the aggregated fuzzy scales a/y (i, j — I, 
2, ... ,n) are transferred into crisp scales a,/ within the range [0, 9]. 



Step 6. Calculate the priority of weights factors 

Let E\, F 2 , . . . , Fn be a set of factors in one section, a,y is the defuzzified scale 
representing the quantified judgment on F, comparing to Fj . Pair-wise comparisons 
between F, and F j in the same section thus yields a n-by-n matrix defined as follows: 





Fi 


r2 


... Fn 


El 


1 


a\2 


■ ■ • (^\n 


A = a/y = F 2 


1 

«12 


1 


■ ■ • Cl2n 


Fn 


1 

a\n 


J_ 

ain 


... 1 



(12.5) 
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where a, 7 = 1, aji — and i, j — 1,2, . . . ,n. 

The priority weights of factors in the matrix A can be calculated by using the 
arithmetic averaging method 



Wi 



1 ” 

1 -^-1 Qij 



( 12 . 6 ) 



where w, is the section weight of f/ . Assume f, has t upper sections at different 
levels in the hierarchy, and is the section weight of the i th upper section 

which contains T, in the hierarchy. The final weight wj of T; can be derived by 



W: 



^FKction 
(=1 



(12.7) 



All individual upper section weights of can also be derived by Eq. 12.6 

to prioritize sections within the corresponding cluster in the hierarchy. 



Step 7. Calculate the fuzzy scores 

When the scores and the priority weights of factors are obtained, the final fuzzy 
scores (FS) can be calculated by 



n 

i=l 



( 12 . 8 ) 



i — 1,2, ..., n 



Step 8. Compare the (FS) values using an outranking method 

Several techniques have been proposed to rank final fuzzy scores [11]. In this case 
the approach proposed by Cheng and Cheng [12] will be used. Let A and B be two 
fuzzy numbers defined as follows: 



fA(x) = 



/f(x), X <niA 
/^( v ), X >ma 



fsix) = 



fg(x), X < mg 
f§(x), X >mB 



(12.9) 

( 12 . 10 ) 



where niA and m g are the mean of A and B. The metric distance between A and B 
can be calculated as follows: 



1 1 

I i9A(y) - 9B(y)fdy + I (g^iy) - g§{y)fdy]-2 



D{A,B) = [ 



0 



0 



( 12 . 11 ) 
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Table 12.1 Criteria and 
subcriteria 



Criteria 


Subcriteria 




cii Efficiency 


Cl Technological 


C12 Reliability 

ci 3 Maturity 

C21 Emissions avoided 


C2 Environmental 


C22 Noise 

C23 Land use 

C31 Investment costs 


C3 Economic 


C32 Operation and Maintenance costs 
C33 Service life 



where g^, g^, g^, and g^ are the inverse functions of f^, fg, and /^ , 
respectively. 

In order to rank fuzzy numbers, let the fuzzy number 5 = 0, then the metric distance 
between A and 0 is calculated as follows: 

1 1 

D(A, 0 ) = [ J(g^(y)fdy + j {gliy)fdy]^ ( 12 . 12 ) 

0 0 

the larger the value of D(A,0), the better the ranking of A. 

According to Cheng and Cheng [12], a trapezoidal fuzzy number A = (ai, a 2 , 
fl 3 , 04 ) can be approximated as a symmetry fuzzy number cr], where /x denotes 
the mean of A, and ct denotes the standard deviation of A. The membership function 
of A is calculated as follows: 



fA(x) = 






ilJ.-\-cr)-x 



if fj. — a < X < fx 
if fx<x<fi + (j 



where /x and cr are calculated as follows: 

2(04 — a\) + 03 — (32 



/X = 



fll + 02 + <33 + tt 4 



(12.13) 



(12.14) 

(12.15) 



if fl 2 = < 23 , then A becomes a triangular fuzzy number, where A = (ai, 02 , 04 ) and 
/X and cr can be calculated as follows: 



(04 — <3l) 



/X = 



“b 2 cz2 4“ ^4 



(12.16) 

(12.17) 
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The inverse functions and , of and respectively are shown as follows: 



SaCv) = + *}' (12.18) 

= (m + ct) - cr (12.19) 



12.1.1 Application 

Let us consider a multi-criteria problem consisting in the selection of a renew- 
able energy project for electric generation. Five alternatives are eligible. Biomass 
(Ai), Geothermal (A 2 ), Hydropower (A 3 ), Solar (A 4 ) and Wind (A 5 ), which will 
be evaluated according to three decision criteria and nine sub-criteria as shown in 
Table 12.1. 

Four decision-makers will evaluate the considered criteria to determine the most 
appropriate renewable energy alternative. Each decision-maker provides a decision 
about his/her judgment as a precise numerical value, a possible range of numerical 
value, a linguistic term or a fuzzy number. A score system is shown in Fig. 12.2. 

Table 12.2 presents the decision-makers’ opinions based on alternatives with 
respect to the related criteria. Then these values are converted into STFNs as defined 
in Eq. 12.1. 

Table 12.3 presents the scores and converted STFN for alternative Ai . The aggre- 
gations of the obtained scores, calculated by Eq. 12.2 for the rest of alternatives are 
shown in Table 12.4. 

The pair-wise comparisons of Technological, Environmental and Economic 
criteria and the corresponding STENs are shown in Tables 12.5, 12.6, 12.7. The 
aggregation of STFNs scales are calculated by Eq. 12.3. 

Then, STEN scale of comparisons are defuzzied by using Eq. 12.4 and by 
using Eq. 12.5, the pair-wise comparisons matrices are obtained. By taking into 
account these pair-wise comparisons matrices and using Eq. 12.6, the weights of the 
sub-criteria are obtained. Table 12.8 shows the the final weights of the criteria cal- 
culated using Eq. 12.7. The final Euzzy scores obtained by using Eq 12.8, are shown 
in Table 12.9. 

In the final step, the fuzzy scores are ranked using the method proposed by Cheng 
and Cheng [12]. The ranking results are shown in Table 12.10. The ranking of 
the alternatives is as follows: A 5 , A 4 , A 2 , Ai and A 3 . According to this ranking, 
alternative A 5 (Wind) the most appropriate alternative. 

12.2 FAHP Chang’s Extent Analysis 

One of the latest approaches on solution processes of EAHP methodology is based 
on Chang’s extent analysis [9, 10]. The Chang’s extent analysis is relatively easy 
while comparing the other approaches on FAHP. The fundamentals of the method 
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Table 12.3 


Scores and converted STFN for Biomass 








DMi 


DM 2 


DM 3 


DM 4 


Cll 


(5 7.5 7.5 10) 


(8 8 8 8) 


(7 7 7 7) 


(8 8 8 8) 


C 12 


(7788) 


(7 7 7 7) 


(8 8 8 8) 


(7 7 8 8) 


Cl3 


(8 8 8 8) 


(8 8 8 8) 


(7 7 7 7) 


(6 7 7 8) 


C 21 


(5 7.5 7.5 10) 


(7 7 7 7) 


(8 8 99) 


(7 8 8 9) 


C 22 


(8 8 8 8) 


(5 7.5 7.5 10) 


(5 7.5 7.5 10) 


(7 7 8 8) 


C23 


(5 7.5 7.5 10) 


(7 7 7 7) 


(8 8 8 8) 


(8 8 8 8) 


C31 


(4 4 4 4) 


(5 7.5 7.5 10) 


(7 7 7 7) 


(8 8 8 8) 


C32 


(5 5 5 5) 


(5 7.5 7.5 10) 


(8 8 8 8) 


(7 7 7 7) 


C 33 


(4 4 4 4) 


(5 7.5 7.5 10) 


(7 7 7 7) 


(7 7 7 7) 


Table 12.4 Aggregated STFNs 


Bioma.ss Geothermal 


Hydropower 


Solar 


Wind 



Cll 


(7 7.6 7.6 8.2) 


(6.5 7.4 7.6 8.5) 


(6.5 7.4 7.6 8.5) 


(6.7 6.7 7.2 7.2) 


(8.5 8.5 8.7 8.7) 


C 12 


(7.2 7.2 7.7 7.7) 


(7.2 7.2 7.7 7.7) 


(7 7 7.2 7.2) 


(7 7 7.2 7.2) 


(3.2 4.1 4.4 5.2) 


C13 


(7.2 7.5 7.5 7.7) 


(8 8 8 8) 


(7.2 7.2 7.5 7.5) 


(7 7 7 7) 


(4.7 5.6 5.9 6.7) 


C 21 


(6.7 7.6 7.9 8.7) 


(6.2 7.1 7.1 8) 


(4.5 4.7 5 5.2) 


(7.2 7.9 8.1 8.7) 


(8.1 8.7 9 9) 


C 22 


(6.2 7.5 7.7 9) 


(7.1 7.7 8 8) 


(6.5 7.1 7.1 7.7) 


(7.7 8.6 8.6 9.5) 


(8.1 8.7 8.7 8.7) 


C23 


(7 7.6 7.6 8.2) 


(6.2 6.9 7.1 7.7) 


(4.7 4.7 5 5) 


(7.5 8.1 8.4 9) 


(8.1 9 9.2 9.5) 


C31 


(6 6.6 6.6 7.2) 


(7 7.9 8.1 9) 


(6.2 7.1 7.1 8) 


(7.5 7.5 7.7 7.7) 


(7.9 8.5 8.7 8.7) 


C 32 


(6.2 6.9 6.9 7.5) 


(6.7 7.4 7.6 8.2) 


(5.2 5.2 5.2 5.2) 


(7.7 7.7 7.7 7.7) 


(7.6 8.2 8.7 8.7) 


C 33 


(5.7 6.4 6.4 7) 


(7.9 8.5 8.5 8.5) 


(6.5 7.1 7.4 8) 


(7.5 7.78 8) 


(8.1 9 9.2 9.5) 



can be examined in four steps: Obtaining priority weights. Comparing degrees of 
possibility, Obtaining the weight vector and Ranking of the alternatives: 

Step 1. Obtaining priority weights 

The first step in this method is to use Fuzzy triangular numbers for pairwise com- 
parison by means of FAFIP scale and next to use extent analysis method to obtain 
priority weights by using synthetic extent values. The value of fuzzy synthetic extent 
with respect to the rth object is represented as: 

m n m 

s = S"i<S>iSS (12.20) 

7=1 /=i 7=1 

and fuzzy addition operation of m extent analysis values can be performed for par- 
ticular matrix such that 

m m m m 

= ( 12 . 21 ) 

,/=l ./=1 ,/=l /=1 

to obtain X7=i ^'gi ■ 
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Table 12.5 Fuzzy aggregation of technological criteria 





Cll 




C\2 




Cl 3 






Scale 


Converted SFTN 


Scale 


Converted SFTN 


Scale 


Converted SFTN 




DMi 




3 a4 


(3 3 4 4) 


3 


(3 3 3 3) 




DM2 




5 


(5 5 5 5) 


2a3 


(2 2 3 3) 




DM2 




5 


(5 5 5 5) 


4 


(4 4 4 4) 




DM4 




4 


(4 4 4 4) 


3a4 


(3 3 4 4) 


Aggregation 


DMi 


(1111) 




(4.3 4.3 4.5 4.5) 


6 


(3 3 3.5 3.5) 
(6 6 6 6) 




DM2 








5 a 6 


(5 5 6 6) 




DM2 








6 


(6 6 6 6) 




DM4 








4a5 


(4 4 5 5) 


Aggregation 


DMi 






(1111) 




(5.3 5.3 5.8 5.8) 



DM2 

DM2 

DM^ 

Aggregation (1111) 



Table 12.6 Fuzzy aggregation of environmental criteria 





C21 




C22 




C 23 






Scale 


converted SFTN 


Scale 


converted SFTN 


Scale 


convertedSFTN 




DMi 




6 a 7 


(6 6 7 7) 


1 a2 


(112 2) 




DM2 




7a8 


(7 7 8 8) 


1 


(1111) 




DM2 




6 a 7 


(6 6 7 7) 


1 a2 


(112 2) 




DM4 




4 


(4 4 4 4) 


3a4 


(3 3 4 4) 


Aggregation 


DMi 


(1111) 




(6.3 6.25 7.25 7.3) 


2 


(1 1 1.75 1.8) 
(2 2 2 2) 




DM2 








2a3 


(2 2 3 3) 




DM2 








1 


(1111) 




DM4 








2a3 


(2 2 3 3) 


Aggregation 


DMi 






(1111) 




(1.8 1.8 2.2 2.3) 




DM2 














DM2 














DM4 












Aggregation 












(1111) 



Then, the fuzzy addition operation of — 1,2, ... ,m) values such that 

n m II fi n 

1 = 1 ,/ = l 1 = 1 /=! i = l 



( 12 . 22 ) 
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Table 12.7 Fuzzy aggregation of economic criteria 





C 31 




C32 




C33 






Scale 


Converted SFTN 


Scale 


Converted SFTN 


Scale 


Converted SFTN 




DMi 




3a4 


(3 3 4 4) 


1 a 2 


(112 2 ) 




DM2 




4 


(4 4 4 4) 


1 


( 1111 ) 




DM2 




3 


(3 3 3 3) 


2 


(2 2 2 2 ) 




DMi 




4 


(4 4 4 4) 


1 a 2 


(112 2 ) 


Aggregation 


DMi 


( 1111 ) 




(3.5 3.5 3.7 3.8) 


5 a 6 


(1.3 1.3 1.7 1.8) 
(5 5 6 6 ) 




DM2 








5 a5 


(5 5 5 5) 




DMi 








5 


(5 5 5 5) 




DMi 








4a5 


(4 4 5 5) 


Aggregation 








( 1111 ) 




(4.8 4.8 5.25 5.3) 



DMi 

DM2 

DMi 

DMi, 

Aggregation (1111) 



Table 12.8 Weight vectors 



Technological 


10.06 


0.03 


0 . 01 ] 


Environmental 


10.34 


0.13 


0.14] 


Economic 


[0.14 


0.10 


0.05] 



Table 12.9 Fuzzy scores 



Table 12.10 Ranking of 
alternatives 



Alternative 


Euzzy Score 


Ai 


(6.5 7.3 7.4 8.2) 


A 2 


(6.7 7.4 7.6 8.2) 


^3 


(5.4 5.8 6 6.3) 


A 4 


(7.4 7.8 8 8.5) 


As 


(7.9 8.5 8.8 8.5) 









Alternative 


D(A,0) 


Ranking 



Ai 


10.47 


4 


^2 


10.52 


3 


^3 


8.34 


5 


A 4 


11.24 


2 


As 


12.04 


1 



are performed to obtain [X/'=i 27= i 

At the end of Step 1, the inverse of the determined vector can be expressed as 
follows: 



1 1 1 . 

zLi/r 



(12.23) 
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Table 12.11 Alternatives and 
criteria 



Alternatives 

A I Wind power 10 < P < 50Mw 
A 2 Hydroelectric 25 < P < 50Mw 
A 3 Solar thermo-electric P > lOMiti 
A 4 Biomass (energetic cultivations) P < 5Mw 
A 5 Biomass (forest and agricultural wastes) P < 5Mw 
A(i Biomass (Co-combustion in conventional central) P > SOMw 
A-] Biogas (electric generation) P < 2Mw 
Criteria 

Power (P){Mw) 

Investment ratio (IR) 

Operating hours (OH) 

Useful life (UL) 

Operation and maintenance costs(OM) 

Emissions avoided (tCo2) 



Table 12.12 EAHP scale 



Definition 


Intensity of importance 


Equal 


(1,U) 


Weak 




Eairly strong 




Very string 


(|>3.2) 


Absolute 


(l4.|) 



Step 2. Comparing degrees of possibility 

The degree of possibility of M 2 = (I 2 , m 2 , M 2 ) > M\ = (/i, mi, mi) is defined as 

y(M 2 > Ml) = sup[min{XjTi(x), X7Z2(y))]y>x (12.24) 

where x and y, are the values on the axis of membership function of each criterion 
and it can be represented as follows: 



V{M2 > Ml) = hgt{Mi n M2) = XTT2{d) = 



1 

0 

(h-u2) 

(m2-U2)-(mi-ul) 



for m 2 > mi 
for /2 > ^1 
otherwise 

(12.25) 



where d is the ordinate of the highest intersection point between firci and ix7t2- 
To be able to compare we need both the values of V (Mi > M 2 ) and V {M 2 > Mi). 



Step 3. Obtaining tbe weight vector 

The degree possibility for a convex fuzzy number to be greater than k convex fuzzy 
numbers M;(/ = 1, 2, A:) can be defined by 
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Table 12.13 Evaluation of RE project choice criteria with respect to goal 



IR 



OH 



UL 



OM 



tco2 



p 


(1,U) 


(f,3,2) 


(2 1 2. 
(i 1 1) 

1 . 7 , 3 , 5 ; 


(i.i.l) 


(2 1 2. 
W’ 3’ 52 




IR 


''7’ 3’ 5'' 
(|> 2 . |) 


(U,l) 


(iil) 


(1,U) 


OH 


(f,3,|) 


(1,U) 


(2 1 2x 

''S’ 2’ 32 




(2 1 2x 

W’ 3’ 52 


UL 


(fA.i) 


(j, 2 , j) 




(1,U) 


(144) 


(U,l) 


OM 


(f,3,7) 


(1,U) 




( 5 . 4 , |) 


(1,U) 


(2 1 2x 

''S’ 2’ 32 


tco2 


(2,2, |) 


(|,1,2) 


(|.3. 2) 


(1,U) 


(2,2.5) 


(U,l) 



V(M > Ml, M 2 , , Mk) =y[(M > Ml), (M > M 2 ), . . . , (M > M*)] 

minV (M > M;) 



(12.26) 



;■ = 1,2, ...,yt 

Assume that, d‘{Ai) = minV {Si > Sk) fork = 1,2, . . . , nv, k ^ i. Then, the 
weight vector is given by 



IT' = {d'{Ax,d‘(A2), d’{An)f (12.27) 

where A/ (/ = 1 , 2, . . . , n) are n elements. Via normalization, the normalized weight 
vector is: 



IT = {d{Ai, d{A2), ..., d{An)f (12.28) 

where IT is a nonfuzzy number that gives priority weights of an attribute or an 
alternative over other. 

Step 4. Ranking of the alternatives 

Based on the weights of the alternatives with respect to the main goal, the ranking 
of the different alternatives is obtained 



12.2.1 Application 

In this subsection the FAHP model is applied to a renewable energy project selection 
for electric power generation. Seven alternatives are available which will be evaluated 
according to six criteria as shown in Table 12.11. The evaluation scale used by experts 
is shown in Table 12. 12. The criteria considered are compared with respect to goal in 
Table 12.13, and comparative evaluations of the alternatives with respect to criteria 
are shown in Tables 12.14, 12.15, 12.16, 12.17, 12.18, 12.19. 



At the end of Step 1, the normalized weight vector obtained from Table 12.13 is 
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Table 12.14 Evaluation of RE project choice criteria with respect to Power 



Ai A2 A3 A4 A5 Ae A7 



Ai 


(1,U) 


(f.i.l) 


( 5 . 1 . 2 ) 


(2.3.2) 


(i.3.2) 


( 3 . 1 . 2 ) 


(u4,|) 


^2 


( 5 . 1 . i) 


(U,l) 


(j.i.i) 


( 2 . 3 . 1) 


(i.3.2) 


(1,1.1) 


(2.4.2) 




(ii^) 


d-ij) 


(l.U) 


( 2 . 3 . 2 ) 


(§.3.2) 


(1,1.1) 


(2.4,2) 


M 


A 1 2^ 

A’ V r 

^7’ 3’ 5 ) 


(i 1 2 x 

A’ V r 

( 7 . 3 . 5 ) 


(2 1 2^ 
A’ V r 

( 7 . 3 , 5 ) 


(l.U) 


(1,1,1) 


a 1 2x 
U. 4’ 1> 

(1 i 

'^9’ 4’ 1> 


(1.1,1) 


As 


(1,1.1) 


(1,1,1) 


(1.1,1) 




dll) 


(l.U) 


(U.l) 


(2,4,2) 


(2.4.1) 


(1,1.1) 


( 2 . 4 , f) 


Ai 


(2 1 2x 

V 9 , 4 , -]) 


(2 1 2^ 
V 9 , 4 , -]) 


(1,1,1) 


(1,1,1) 


/2 1 2x 

4’ 1’ 


(1,1,1) 



Table 12.15 Evaluation of RE project choice criteria with respect to Investment Ratio 





Ai 


A 2 


A 3 


A4 


As 


A6 


A? 


Ai 


(1,1,1) 


(2 1 2. 
'•S’ 2 ’ 3^ 


(2.4.?) 


(f.3.2) 


(f,3,2) 


(1,1,1) 


(1,1,1) 


A 2 


(2.2.|) 


(1,1,1) 


(2.4.!) 


(2.3,2) 


(i.3.2) 


(2.2. i) 


( 2 . 2, 2 ) 


A3 


( 5 . 4 . 7 ) 


a 1 2^ 

'• 9 ’ 4 ’ 1’ 


(1,1,1) 


( 7 . 5 . 5 ) 


( 7 . 3 . 5 ) 


( 5 .U 7 ) 


(i.2,f) 


A4 


3’ 5’ 


'•7’ 3’ S'' 


(1.3. 2) 


(1,1,1) 


(1,1,1) 


W’ 3’ 5'' 


(1.1,1) 


As 


(2 1 2x 

3’ 5’ 


(2 1 2^ 
'•7’ 3’ 5'' 


(1.3.2) 


(1,1,1) 


(1.1,1) 


(2 1 1) 
''7’ 3’ S'' 


(1,1,1) 


^6 


(1.1,1) 


U, 2’ 3'' 

(2 1 U 

'•S’ 2’ S'' 


(2.4.^) 


(2.3,2) 


(2.3.2) 


(1,1,1) 


(1,1,1) 


A7 


(1,1,1) 


(iii) 


(1,1,1) 


(1.1,1) 


(1,1,1) 


(1,1,1) 


Table 12.16 Evaluation of RE project choice criteria with respect to Operating Hours 




Ai 


^2 


A3 


A4 


As 


A6 


Ay 



Ai 


(1,1,1) 


(i.2.i) 


.2 1 2x 

2’ S-' 


(2 1 2. 
U, 4, y 7 


(2 1 2. 

Vo, 4, j ) 

(1 1 2x 

V9, 4, y 7 


(2 1 2^ 
Vo, 4, j ) 

(1 i i ) 

V9, 4, j ) 


(5. J.7) 


A2 


(5.2.3) 


(1,1,1) 


'^S’ 2’ 3 '' 


'• 9 ’ 4 ’ 7 '' 


(9. 4, 7) 


A3 


(i. 2 . 5 ) 


(i.2.f) 


(1,1,1) 


(2 1 2^ 
'• 9 ’ 4 ’ 7 '' 


.2 1 2^ 
V9, 4, j) 


a 1 2x 
V9, 4, 7; 


a 1 2x 
'• 7 ’ 3 ’ S'' 


A4 


(2.4,^) 


( 2 , 4 . 2 ) 


( 2 . 4 , 2 ) 


(1,1,1) 


(1,1,1) 


(1,1,1) 


(1,1,1) 


As 


(2.4,^) 


( 2 , 4 . 2 ) 


( 2 . 4 , 2 ) 


(1,1,1) 


(1,1,1) 


(1,1,1) 


(1,1,1) 


Ae 


(2.4,^) 


(^. 4 .^) 


(2.4,^) 


(1,1,1) 


(1,1,1) 


(1,1,1) 


(1,1,1) 


A4 


( 2 . 4 , 2 ) 


( 2 . 4 . 2 ) 


(f. 3 , 2 ) 


(1,1,1) 


(1,1,1) 


(1,1,1) 


(1,1,1) 



Table 12.17 Evaluation of RE project choice criteria with respect to Useful Life 



Ai A2 A3 



Ai 


(1,1,1) 


(1 1 2 x 
'•S’ 2 ’ 3 '' 


'•S’ 2 ’ 


2) 

3 '' 


A2 


(i. 2 ,i) 


(1,1,1) 


(1,1,1) 




A3 


(2.2,2) 


(1,1,1) 


(1,1,1) 




A4 




a 1 2 x 
'• 7 ’ 3 ’ S'' 

(i 1 2 x 

'• 7 ’ 3 ’ S'' 


'• 7 ’ 3 ’ 


|) 


As 


(5. 2. 5) 


'• 7 ’ 3 ’ 


5) 


Ae 


(1,1,1) 


/2 1 2 x 

'•S’ 2 ’ 3 '' 


V- - 

*• 3 ’ 2 ’ 


3) 


Ay 


(1,1,1) 


(2 i 2) 
15.2.31 


(2 2 
I5. 2’ 


1 ) 



A4 


As 


Ae 


Ay 


(2.2. i) 


(2.2. i) 


(1,1,1) 


(1,1,1) 


(f, 3 . 2 ) 


(2.3. 1 ) 


(i. 2 ,i) 


(2. 2, 2) 


(2.3.2) 


(|. 3 . 2 ) 


(i. 2 . 2 ) 


(7. 2, |) 


( 1 , 1 , 1 ) 


( 1 , 1 , 1 ) 


(f.i,|) 


(f,i,|) 


( 1 , 1 , 1 ) 


( 1 , 1 . 1 ) 


(5. 2 . 5) 


(5. 2. 3) 


(2. 2, 2) 


(2.2. f) 


(1,1,1) 


(1,1,1) 


(7. 2, 7) 


(^. 2 . 5 ) 


(1,1,1) 


(1,1,1) 
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Table 12.18 Evaluation of RE project choice criteria with respect to Operation and Maintenance 
costs 





Ai 


42 


43 


A4 


45 


46 


47 


Al 


(U.l) 


(2 1 2. 
'•S’ 2’ 3^ 


(i.3.2) 


(2 1 2. 
''S’ 2’ 3^ 


(2 1 2. 

''S’ 2’ 3^ 


(2 1 2^ 
'•S’ 2’ 3^ 


(1,1,1) 


A2 


(i.2.5) 


(1,U) 


(2.4.2) 


(1,U) 


(l.U) 


(1,1,1) 


(i.3.2) 


A 3 


ty, 3, 


(2 1 2^ 
'•9’ 4’ 


(U.l) 


(2 1 2x 

''9’ 4’ 7’ 


(2 1 2^ 
''9’ 4’ 7^ 


(2 1 2^ 
'•9’ 4’ 7> 


(5.4.7) 


A4 


(|.2. |) 


(l,U) 


(2.4.2) 


(1,U) 


(l.U) 


(1,1,1) 


(i.3.2) 


45 




(1,U) 


( 5 . 4 .^) 


(1,U) 


(l,U) 


(1,1,1) 


(|,3,2) 


^6 


{j, 2, j) 


(1,U) 


a. 4.!) 


(U.l) 


(U.l) 


(1,1,1) 


(i.3,i) 


Ay 


(U,l) 


(2 1 2) 
3’ 5> 


(f,3. 2) 


(2 1 2x 

3’ 52 


(2 1 2) 
'7' 3’ 5I 


(2 1 2, 
I7. 3. 5I 


(1,1,1) 



Table 12.19 Evaluation of RE project choice criteria with respect to Emissions avoided 





4 i 


2^2 


2^3 


A4 


45 


2^6 


Ay 


4 i 


(1,1,1) 


(5.4, |) 


(7.3.^) 


(3.3.^) 


(i. 3 . 3 ) 


(3.2,?) 


(3.4, |) 


42 


(2 1 2v 

'• 9 ’ 4 ’ 1 ’ 


(1,1,1) 


(5.3.3) 


(f if) 


(7.^5) 


(5.4.7) 


(1,1,1) 


2^3 


(1 1 2>j 

'• 7 ’ 3 ’ 5 ^ 


(i. 2 .i) 


(1,1,1) 


(2 1 2. 
W’ 3 ’ S-^ 


(2 1 
'• 7 ’ 3 ’ S2* 


(2 1 2. 
'• 9 ’ 4 ’ 7 ’ 


( 7 . 2 ,i) 


A4 


(7.3.5) 


(i. 3 . 7 ) 


(i, 3 . 3 ) 


(1,1,1) 


(1,1,1) 


/2 1 2>j 

'•S’ 2 ’ 32* 


(i. 3 . 3 ) 


45 


(5. 7. 3) 


(?. 3 , 2 ) 


(?, 3 . 3 ) 


(1,1,1) 


(1,1,1) 


(2 1 2-j 

'•S’ 2 ’ 3 '' 


(?, 3 , 3 ) 


2^6 


isJ) 


(5.4, i) 


(3.4.!) 


( 5 . 3 . 3 ) 


( 5 . 3 . 3 ) 


(1,1,1) 


(§. 3 . 3 ) 


Ay 


(- - -) 
'• 9 ’ 4 ’ 7 ’ 


(1,1,1) 


(5.1-3) 


(2 1 2^ 
I7. 3. 5I 


(2 1 2x 

I7. 3. 5I 


(2 1 2^ 
I7. 3. 5I 


(1,1,1) 



iZ2Xr' = ( 

i=\ j=\ 



1 



1 



1 



55.22 46.08 38.37 



) 



In order to identify the computation stages clearly, the pair-wise comparisons 
from Table 12.13 are evaluated as follows: 



Sp = (5.25, 6.33, 7.73) (g) ( 



SjR = (3.64,4.17,4.97) ® ( 



SoH = (6.35, 7.83, 9.57) g) ( 



Sul = (5.89, 7.25, 8.79) g( 



SoM = (9.07, 10.5, 12.17) g ( 



Steal = (8.17, 10.0, 12.0) g ( 



1 



1 



1 



55.22 46.08 38.37 
1 1 1 

55.22 46.08 38.37 
1 1 1 

55.22 46.08 38.37 
1 1 1 

55.22 46.08 38.37 
1 1 1 

55.22 46.08 38.37 
1 1 1 

55.22 46.08 38.37 



) = (0.10,0.14,0.20) 
) = (0.07,0.09,0.13) 
) = (0.12,0.17,0.25) 
) = (0.11,0.16,0.23) 
) = (0.16,0.23,0.32) 
) = (0.15,0.22,0.31) 



After determining these results the following are obtained 
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Table 12.20 Weights of alternatives with respect to goal 





p 

(0.09) 


IR 

(0) 


OH 

(0.18) 


UL 

(0.15) 


OM 

(0.30) 


tco2 

(0.28) 


Ai 


0.22 


0.27 


0 


0.07 


0 


0.39 


42 


0.22 


0.46 


0 


0.40 


0.25 


0 


^3 


0.22 


0 


0 


0.40 


0 


0 




0 


0 


0.25 


0 


0.25 


0.18 


45 


0 


0 


0.25 


0.07 


0.25 


0.26 




0.33 


0.27 


0.25 


0.07 


0.25 


0.26 


Ay 


0 


0 


0.25 


0.07 


0 


0 



V{Sp > Sir) = 1 
V(Sp > Soh) = 0.72 
ViSp > Sul) = 0.81 
V (Sp > Sqm) = 0.31 
V(Sp > Sre) = 0.38 

Applying Eq. 12.26: 



minC^p > Sir, Sjp, Soh, Sul, Sqm) — 0.31 
miniSiR > Sp, Sjp, Soh, Sul, Som) — 0 
minC^/p > Sp, Sir, Soh, Sul, Som) — 0.60 
> Sp, Sir, Sip, Sul, Som) = 0.50 
minC^f/L > Sp, Sir, Sip, Sqh, Som) — 1 
min(5'oM > Sp, Sir, Sip, Son, Sul) — 0.94 

Therefore, the weight vector from Table 12.13is W‘ — (0.31, 0, 0.60, 0.50, 1, 0.94)^ 
and the normalized weight vector is VT(0.09, 0, 0.18, 0.15, 0.30, 0.28)^. The same 
systematic approach is considered for other evaluations and the weights obtained 
from Tables 12.14, 12.15 12.16, 12.17, 12.18, 12.19 are: 

Table 12.4 (0.22,0.22, 0.22,0, 0,0.33,0)^ 

Table 12.5 (0.27, 0.46, 0, 0, 0, 0.27, 0)^ 

Table 12.6 (0, 0, 0, 0.25, 0.25, 0.25, 0.25)^ 

Table 12.7 (0.07, 0.4, 0.4, 0, 0, 0.07, 0.07)^ 

Table 12.8 (0, 0.25, 0, 0.25, 0.25, 0.25, 0.25)^ 

Table 12.9 (0.39, 0, 0, 0.18, 0.18, 0.26, 0)^ 



Finally, weights of the alternatives with respect to the goal are given in Table 
12.20 and ranking of the alternatives based on overall results is given in Table 12.21. 
As the result of the evaluation, the three Biomass alternatives are the best ranked by 
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Table 12.21 Ranking of 
alternatives 



Alternative 



^6 


0.23 


A 4 


0.17 


A 5 


0.17 


^2 


0.15 


hi 


0.14 


^3 


0.08 


Al 


0.06 



the method, and the Co-combustion in conventional central plant (P > 50Mw) is 
selected as the best one with a priority weight of 0.23. 



12.3 Concluding Remarks 

When using FAHP the decision-maker is only asked to give judgments about either 
its preference of one alternative on one criterion against another or the relative impor- 
tance of one criterion against another. However, when the number of alternatives and 
criteria grows, the method requires complex computations and the process becomes 
cumbersome. The method has been used by Wang et al. [13] for trigeneration systems 
selection and evaluation and by Jaber et al. [14] to evaluate space heating systems 
used in Jordan based on the fuzzy sets and analytical hierarchy process. Kahraman 
et al. [15] use two fuzzy multicriteria decision-making methodologies, FAHP and 
Fuzzy axiomatic design, to select the most appropriate renewable energy alternative 
for Turkey, Kaya and Kahraman [16] propose and integrated VIKOR-FAHP in order 
to determine the best renewable energy alternative for Istanbul. Talinli et al. [17] 
present a comparative analysis of three different energy production process (EPP) 
scenarios for Turkey using FAHP. 
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Chapter 13 

The Shapley Value 



In fuzzy Multi-criteria decision-making problems, the ranking of alternatives must 
take into account their fuzzy scores in all criteria, the weights assigned to each 
decision criterion, the possible difficulties of comparing two alternatives when one is 
significantly better than the other on a subset of criteria, but much worse on at least one 
criterion from the complementary subset of criteria and the decision-maker’s attitude 
toward the risk associated with evaluation. Therefore, the relationships among criteria 
are crucial for adequate treatment of fuzzy decision-making because they refiect the 
structure of interaction among the criteria and represent decision-maker’s preferences 
of the criteria. The method presented in this chapter produces an evaluation of the 
marginal contribution of each of the decision criteria to the overall goal refiecting 
the relationships among the decision criteria and decision-makers’ preferences. 



13.1 The Shapley Value 

The global importance of a particular criterion is not solely determined by that par- 
ticular criterion but also by the value of all other criteria considered in the evaluation 
process [1]. Let us consider a set of decision criteria as C = [ci, C 2 , . . . , c„, ] and /x 
a fuzzy measure on C such that 

p,(C) = ^/xKq) for / = 1,2, ...,n (13.1) 

i 

with /x(<t>) = 0 and /x(C) = 1, where <I> = null set and /x(C) is the weight or the 
importance value of the criterion c, . To determine the expected marginal contribution 
of a particular criterion (player) to the overall goal, the importance index or Shapley 
value of criterion c,- with respect to /x is defined as: 

Z (Af- A)!(A- 1)! 

^ ^[/x(A) - /x(A - Cl)] (13.2) 

where N is the number of decision criteria and A is any combination of decision 
criteria containing criterion c, . Thus, in addition to the usual weights on criteria 
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Table 13.1 Fuzzy numbers 
for linguistic variables 



Definition 


Fuzzy number 


VG (Very Good) 


(0.8 0.9 I.O I.O) 


G (Good) 


(0.6 0.7 0.8 0.9) 


AA (Above Average) 


(0.5 0.6 0.7 0.8) 


A (Average) 


(0.4 0.5 0.5 0.6) 


BA (Below Average) 


(0.2 0.3 0.4 0.5) 


P/LI (Poor/Low Important) (0.0 0.0 0. 1 0.2) 



taken separately, weights on any combination of criteria must also be debned in Multi- 
criteria Decision-Making problems. The method proposed by Shapley [2] consists 
of the following steps. 

Step 1. Establishing weights for subcriteria 

Consider a decision problem with criteria and subcriteria. First, the decision-maker 
expresses, as fuzzy numbers, the importance weights assigned to each of subcriteria 
according to Table 13.1. Next, these values are defuzzied and then normalized. 

Step 2. Calculating the Shapley value for subcriteria 

In order to calculate the Shapley value for subcriteria (Eq. 13.2), decision-makers are 
asked to assign the weight to the combinations of subcriteria, based on the weights 
established for subcriteria. Thus, w{c\\, C 12 ) represents the average weight assigned 
to the combination of subcriteria 1 and 2 and w(<l>) is the weight to the null set 
of subcriteria. These subcriteria do not compensate each other, that is, a bad score 
on one criterion cannot be adjusted by a good score on any other complementary 
criteria. However, in practice, these subcriteria by and large corporate each other, 
that is, an increase (or decrease) in the degree to which one criterion is satished often 
increases (or decreases) the degree to which another criterion is satished. Therefore, 
the weight assigned to the combination of any two criteria must be at least equal to 
the sum of their individual weight assigned separately [2]. That is, 

w(c]_, C2) > w(ci) + w(c2) (13.3) 

Step 3. Rating of performance of alternatives 

To know the scores for different alternatives against each of the different criteria, 
alternative’s performances on each of the subcriteria are to be rated by the decision- 
makers. 

Step 4. Calculating the Shapley value for decision criteria 

Decision-makers are asked to assign the weights to each of the decision criteria and 
combinations of them, then the Shapley value of each decision criterion is calculated 
using Eq. 13.2. The values indicate the expected marginal contribution of global 
importance of each decision criterion to the overall goal. 
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Table 13.2 Criteria and 
subcriteria 



Table 13.3 Ranking of 
subcriteria 



Criteria 


Subcriteria 


Cl Economic 


cii Investment cost 
C12 Net present value 


C2 Technological 


Ci3 Payback period 
C21 Reliability 
C22 Efficiency 


C3 Environmental 


C23 Maturity 

C31 Emissions avoided 

C32 Noise 




C33 Land use 



Economic 


Environmental 


Technological 


cii VG 


^'21 


G 


C31 G 


Cl2 G 


t'22 


G 


C32 A 


Cl3 G 


C23 


AA 


C33 G 



Step 5. Ranking of alternatives 

In order to rank the alternatives, the Total Score (TS) for an alternative can be cal- 
culated as follows [3]: 



75 = ^;ri../x(c*) fork = 1,2, :,n (13.4) 

where /x(ci) is the Shapley value on subcriterion k, and is the score of the alter- 
native against the subcriterion k. 



13.2 Application 

In this section, in order to show the use of the Shapley value, the selection of a 
renewable energy project for electric generation is considered. Four alternatives are 
eligible: Solar (Ai), Wind (A 2 ), Hydropower (A 3 ) and Biomass (A 4 ). The four 
alternatives will be evaluated taking into account three decision criteria (Economic, 
Technological, and Environmental) which consists of the subcriteria shown in Table 
13.2. The importance weights assigned by the decision-maker to each of these sub- 
criteria are shown in Table 13.3. 

In order to calculate the crisp scores (defuzzied value of a trapezoidal fuzzy 
number) for subcriteria, the following equation was proposed by Kaufmann and 
Gupta [4]: 



(xi -f X2 -F X3 -t- X4) 



e — 



4 



(13.5) 
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Table 13.4 Crisp scores for subcriteria 



Economic 




Technological 




Environmental 




Subcriterion 


Crisp score 


Subcriterion 


Crisp score 


Subcriterion 


Crisp score 


Cll 


0.925 


C2\ 


0.750 


C31 


0.750 


C12 


0.750 


C22 


0.750 


C32 


0.500 


C13 


0.750 


C23 


0.650 


C33 


0.750 



Table 13.5 Normalized 
weights for subcriteria 



Economic Technological Environmental 

Subcriterion weight Subcriterion weight Subcriterion weight 



Cll 


0.381 


C 21 


0.349 


C 31 


0.375 


Cl2 


0.309 


C 22 


0.349 


C 32 


0.250 


Cl3 


0.309 


C 23 


0.302 


C 33 


0.375 



Table 13.6 Weights assigned 
to different combinations of 
subcriteria 



Economic 




Technological 


Environmental 


Combination weight Combination weight Combination weight 


>C(C11,C12) 


0.80 


>C(C21,C22) 


0.75 


>c(C3l,C32) 


0.80 


u;(cii, C13) 


0.75 


IC(C21. C23) 


0.80 


IC(C31, C33) 


0.75 


IC(C12, C13) 


0.80 


W(C 22 , C23) 


0.75 


IC(C32, C33) 


0.80 



The crisp scores (defuzzied valued) for subcriteria are shown in Tables 13.4 and 
13.5 shows the normalized weights for each subcriterion obtained by dividing scores 
of each dk by XLi Cii- 

Let us assume that the weights assigned to different combinations of subcriteria 
are as presented in Table 13.6. 

Using Eq. 13.2 the Shapley value for subcriteria can be calculated as follows: 

1. When A is cn alone 



(3- 1)!(1 - 1)! 

m(cii)i = ^ ^[w(cii) - «;(<!>)] = 0.127 

2. When A is the combination of c\ i and cu 

(3-2)!(2- 1)! 

/x(cii )2 = [w(cii, C 12 ) - w(ci, 2 )] = 0.082 

3. When A is the combination of c\\ and C 13 

(3-2)!(2- 1)! 

/x(cii )3 = g, [w(cii, C 13 ) - w(ci, 3 )] = 0.073 

4. When A is the combination of cu , C 12 and C 13 



(3 - 3)!(3 - 1)! 

/x(cn)4 = [w(cii, C12, C13), -w(ci,2, C13)] = 0.066 



3 ! 
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Table 13.7 Shapley value for 
subcriteria 



M(cii) = 0.349 m(c2i) = 0.349 = 0.346 

/x(ci2) = 0.338 m(c 22) = 0.324 /r(c32) = 0.308 

M(ci3) = 0.313 m(c23) = 0.326 m('^33) = 0.346 



Table 13.8 Alternative’s performance for subcriteria 





Subcriterion 






Subcriterion 




Subcriterion 




Cll 


Cl2 


Cl3 


C21 


C22 


C23 


i?31 


C32 


C33 




VG 


G 


G 


G 


A 


A 


AA 


G 


G 




AA 


VG 


AA 


G 


AA 


VG 


G 


AA 


VG 


^3 


G 


G 


A 


AA 


VG 


AA 


G 


A 


G 


Aa 


VG 


G 


AA 


G 


VG 


G 


G 


VG 


AA 



Table 13.9 Crisp scores 




Economic 




Technolo 


gical 




Environmental 






Cll 


Cl2 


C13 


C21 


C22 


C23 


C31 


C32 


C33 


Ai 


0.925 


0.750 


0.750 


0.750 


0.500 


0.500 


0.650 


0.750 


0.750 


^2 


0.650 


0.925 


0.650 


0.750 


0.650 


0.925 


0.750 


0.650 


0.925 


As 


0.750 


0.750 


0.500 


0.650 


0.925 


0.650 


0.750 


0.500 


0.750 


Aa 


0.925 


0.750 


0.650 


0.750 


0.925 


0.750 


0.750 


0.925 


0.650 



Now, Shapley value for subcriterion c\\ can be obtained as: 

/x(cii) = /x(cii)i + /x(cii )2 + /r(cii)3 + /x(cii)4 = 0.349 

Similarly, Table 13.7 shows the Shapley value for the rest of subcriteria 

Table 13.8 shows the linguistic assessment of the alternative’s performance on 
each of the subcriteria and Table 13.9 shows the crisp scores on each of the subcriteria 
obtained using Eq. 13.3. 

Total scores for alternatives on criterion ci are calculated as follows: 

"0.925,0.750,0.750 
0.650, 0.925, 0.650 
0.750, 0.750, 0.500 
0.925,0.750,0.650 

Thus, Table 13.10 shows the scores for different alternatives against each of three 
different criteria. 

Let us assume that the weight assigned by the decision-maker to each of the 
decision criteria and combinations of them for this evaluation process are the weights 
shown in Table 13.11. The Shapley value of each decision criterion calculated using 
Eq. 13.2 is shown in Table 13.12. 

Scaling these values by a factor N = 3, an importance index greater than 1 
indicates a criterion more important than the the average. From Table 1 3. 12, it is clear 



[0.349,0.338,0.313] = [0.811,0.743,0.672,0.780] 
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Table 13.10 Score matrix for 
alternatives 



Alternatives 


Cl 


C 2 


C 3 


Al 


0.811 


0.587 


0.715 


A 2 


0.743 


0.775 


0.780 


^3 


0.672 


0.739 


0.673 


A 4 


0.780 


0.807 


0.769 



Table 13.11 Weights 
assigned to decision criteria 
and combinations of them 



u;(ci)=0.4 u)(ci, C 2 ) = 0.70 u;(<I> = 0) 

u;(c 2 ) = 0.3 u)(ci, C3) = 0.85 u;(ci, C2, C3) = 1.0 

u;(c 3 ) = 0.3 u)(c 2 , C3) = 0.90 



Table 13.12 Shapley value 
for decision criteria 



Shapley value 0.325 0.300 0.375 

Scaled shapley value 0.975 0.900 1.125 



Table 13.13 Ranking of 
alternatives 



Score 0.783 0.766 0.708 0.692 

Rank A 4 A 2 A\ A 3 



that the decision-maker attaches more weight to criterion 3 than any other criterion 
in the selection process. Using Eq. 13.4, the overall score for each alternative can he 
calculated 

"0.811,0.587,0.715 
0.743,0.775,0.780 
0.672,0.739,0.673 
0.715,0.807,0.769 

Hence, the overall performance score for alternatives and their ranking are as 
shown in Table 13.13. Alternative A4, the Biomass plant is the best ranked and 
therefore, the most preferred alternative. 



[0.325,0.300,0.375] = [0.7080.7660.6920.784] 



13.3 Concluding Remarks 

Historically, the first application of the Shapley value was cost sharing or budgeting in 
situations in which multiple users contribute resources and share the costs or profits of 
the ensuing activities, but since then the method has been used by researchers in cost 
allocation problems [5], water resource allocation [6] and to evaluate the efficiency 
in the exploitation of hydroelectric plants [7]. In energy and environmental analysis, 
the method has been used by Albretch et al. [8], Ang et al. [9] and Wu et al. [10]. 
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